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ABSTRACT 


This final report documents the results achieved during a study of 
the Detail Content of Apollo Orbital Photography under contract NASW-2434 
with NASA Headquarters. This study spanned an 18 month period and was composed 
of three tasks whose objectives were to assess the effect of anomalous V/H 
(Velocity to Height ratio) sensor operation or photoreproduction processes 
upon the detail content of lunar surface photography obtained from the orbiting 
Command Module. This report includes data and conclusions obtaiend from the 
Apollo 15, 16 and 17 missions. 

The specific tasks undertaken included (1) an evaluation of the 
residual motion smear in Apollo 15 panoramic photography due to intermittent 
V/H sensor operation and (2) intercomparison of the detail content of multiple 
generations of Apollo 16 and 17 photography. 

During the Apollo 15 mission the panoramic camera V/H sensor did not 
operate properly and resulted in anomalously high or pre-set nominal command 
V/H rates. In some cases a transition between the anomalous rate and the 
nominal rate occurred during the exposure of a frame. Two such frames (316 
and 360) were examined to deteimiine the subsequent effect on the detail 
content. The evaluation showed that a loss in resolution due to image motion 
blur and exposure differences within a frame could occur due to the inter- 
mittent V/H sensor operation. The lack of correct compensation reduced the 
resolution by 1/3 (to about 6 meters on the lunar surface). Only when the 
commanded V/H is equal to the nominal value can the full resolution potential 
of the panoramic camera be realized; but this may not always occur. Consequently 
it is concluded that a significant portion of the Apollo 15 panoramic imagery 
may have lower resolution than the nominal 2 meter value. The user of this 
imagery should refer to the telemetry data to determine the V/H command 
situation for the frames he is exploiting and account for the effects of 
poor resolution if necessary. 
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An assessment of the detail content of multiple generations of 
Apollo 16 and 17 photography was made to determine the degradation, if any, 
in the quality of the photographic products furnished to various users. 

The original flight film is copied to produce a Master Positive (2P) or 
Direct Negative (2N) copy. These are the earliest generations which the 
user may exploit and only a limited number of such copies are produced. Most 
investigators are likely to use 3rd Generation or 4th Generation Positive 
copies. Unfortunately 3P and 4P copies of the same type of photography were 
not available thus preventing a direct comparison of the most likely user 
products. The evaluation showed that the tone quality of the Master Positive 
is superior to that of the Direct Negative for the panoramic photography. 

This superiority results from a greater dynamic range which means that the 
2P photography contains higher contrast images over a wider range of average 
exposure levels than the 2N photography. The data indicates that this 
difference is preserved in the higher generation products which the user 
is likely to employ. Tonal quality differences were also present in 2nd 
Generation copies of metric photography, however, it is not clear that these 
differences are preserved in the user products. MTF measurements were made 
in order to compare the fine detail content of the various generations of 
photography evaluated. Approximately half of these frames showed some loss 
in fine detail content which is attributed to the reproduction process. An 
increase in the ground resolution element by a factor of 1.5 was observed; 

3 meters for the panoramic photography and 30 meters for the metric photography. 

-This study successfully demonstrated the development and application 
of image evaluation methods for orbital photography of the lunar surface. 
Descriptions of the various elements of software used in this study are 
included in this report . It is recommended that the development of these 
techniques be continued for application to future manned and unmanned space- 
craft involved in planetary exploration. 
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1 . INTRODUCTION 


This final report documents the results achieved during a study of 
the Detail Content of Apollo Orbital Photography under Contract NASW-2434 
with NASA Headquarters. This study was composed of three tasks whose 
objectives were to assess the effect of a V/H anomaly or photo reproduction 
process upon the detail content of lunar surface photography obtained from 
the orbiting Command Module. The report includes results from the photo 
evaluation of the Apollo 15, 16 and 17 missions. An earlier report presented 
the results of a similar study of photography from the Apollo 8, 12, 14 and 
15 missions. This report also contains a detailed description of the Edge 
Gradient Spectral Analyses (EGSA) computer program that was employed 
extensively during these studies. 

The specific tasks completed include (1) an evaluation of the residual 
motion smear in several Apollo 15 panoramic photographs due to intermittent 
V/H sensor operation and (2) intercomparison of, the detail content of multiple 
generations of Apollo 16 and 17 photography. . 

Before proceeding we will briefly review the photography obtained 
during these Apollo missions and its reproduction sequence. In all these 
missions the Scientific Instrument Module (SIM) Bay contained two cameras; 
an Optical Bar Panoramic Camera and a Mapping (or Metric) Camera. The 
Panoramic Camera had a 610 mm (24 inch) focal length lens and employed EK 
Type 3414 film yielding an approximate resolution capability of 130 cycles/mm. 
The Mapping Camera had a 76 mm (3 inch) focal length lens and employed EK 
Type 3400 film. Its resolution capability was about 90 cycles/mm. Both 
cameras contained internal image motion compensation (IMC) mechanisms. 

The effect of erratic operation of the Apollo 15 panoramic IMC was assessed 
in the study. 

The original flight film is reproduced for distribution to potential 
users. Figure 1 shows the reproduction, sequence and terminology used. The 
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flight film is copied to produce a Second Generation or Master Positive (2P) 
copy. Only a limited number of such copies are generated. They are sub- 
sequently copied twice to produce a: fourth Generation (4P) copy which many 
investigators are likely to use. The flight film is also copied on to reversal 
film to obtain a Second Generation Direct (2N) Negative. This is copied once to 
produce a Third Generation (3P) Positive which is also available to investi- 
gators. This sequence eliminates one copy step in the hope of improving 
product quality and reducing production time. An assessment was made of the 
effect of both reproduction sequences on the detail content of Apollo 16 
and 17 photography. 



Figure 1 APOLLO ORBITAL PHOTOGRAPHY REPRODUCTION SEQUENCE 
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2. EVALUATION TECHNIQUES 

All of the methods used in this study to evaluate the Apollo orbital 
photography were developed during previous contracted studies for NASA and 
others. These methods are briefly described in this section of the report. 

In some cases complex computer programs are required. Their operation is 
discussed and the .software detailed in the appendix. 

T^e evaluation of the quality of Apollo Orbital photography can 
require the measurement of various properties of the imagery. The properties 
which define image quality can be divided into four general categories: 

: (1) "resolution" or fine detail content, (2) contrast or tone quality, (3) 
noise level and (4) metric quality. Overall performance can be based upon 
a composite measure which- includes several (or all) of the categories but 
must be determined with the intended use of the imagery in mind. 

The modulation transfer function (MTF) is often used as a measure 
of image resolution or fine detail content. In this study we used this 
image quality measure to determine the image blur resulting from erratic 
operation of the Apollo 15 panoramic camera IMC and to evaluate the nominal 
resolution in multiple generations of Apollo 16 and 17 metric and panoramic 
photography. 

Contrast or tonal quality can be described as the fluctuation in 
density due to large area intensity differences in the scene (i.e. low 
frequency brightness variations) . Such fluctuations are frequently used to 
classify areas or delineate regional boundaries. The most common measure 
of tonal quality is the Hurter-Driffield or D-Log E response curve. It 
describes how the apparent scene brightness is reproduced as density in the 
photograph. The D-Log E response curve was measured for each generation of 
Apollo 16 and 17 imagery evaluated to assess the effect of the reproduction 
sequence on tonal quality or large scale detail content. 
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2.1 Frame Selection - In selecting the frames to be analyzed the 
primary consideration is the availability of test targets for estimating the 
MTF. A shadow-to-sunlight edge inside craters will occur in the low sun 
angle photography. If the sun angle is too high this edge either does not 
exist or is too close to the rim of the crater to be useful. Similarly 
at extremely low sun angles the edge will be sufficiently close to the far 
rim of the crater to make edge analyses impractical. Sun elevation angles 
between 10° to 15° yield the most satisfactory shadow-to-sunlight edge. 

Table 1 indicates the frames analyzed during the study. 

Table 1 

APOLLO FRAMES SELECTED AND ANALYIZED 


MISSION 

NO. 

CAMERA 

FRAME 

NO. 

GENERATIONS 

ANALYIZED 

OBJECTIVE 

15 

PANORAMIC 

316 

360 

4P 

4P 

V/H ANOMALY 

16 



2P, 2N, 3P 
2P, 2N, 3P 

REPRODUCTION 

EVAULATION 

17 

PANORAMIC 

1633 

3111 

2P, 2N, 3P 
2P, 2N, 3P 

MAPPING 

180 

2923 

2P, 2N,4P 
2P, 2N,4P 


2.2 Calculating the NfTF from. Edges . It is common practice to use 
edges to calculate the MTF of an optical imaging system both in the laboratory 
and operationally. The MTF can be calculated using any target by dividing 
the measured image spectrum by the known spectrum of the object. In the 
case of an edge gradient analysis no correction for the target is necessary. 
The derivative of a sharp edge (i.e. gradient) is a narrow pulse (Dirac delta 
function) whose Fourier transform or spectrum is a constant proportional 
to the peak-to-peak amplitude of the edge. In selecting the edge for analysis 
we only require that it be much sharper than the degradation introduced by the 
optical system being evaluated. In the present case the shadow-to-sunlight 
edge sharpness is determined by the penumbra effect produced by the finite 
angular size of the sun. This effect is proportional to the diameter of the 
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crater. To insure sufficient sharpness we selected craters about 10 
resolution elements wide (about 80 microns on panoramic imagery) . 

The edges are scanned using a microdensitometer and recorded on a 
strip chart or magnetic tape. The edge sharpness requirement produces 
targets that are short and curved. Consequently the microdensitometer scanning 
aperture must be small in both dimensions. If the edges were long, a long 
but narrow slit could be used to minimize the noise in the edge trace. Since 
this is not possible the density traces on the strip chart were hand-smoothed 
and digitized using a Calraa digitizer. The data flow is summarized in Figure 2. 
A special software package, not; shown in the figure, supports the Calma to 
generate a card deck containing the edge traces. The Edge Gradient Spectral 
Analyses (EGSA) software package shown in the figure consists of two major 
programs. The initial program reads the digitized density trace, converts 
it to exposure (illuminance trace) and writes an output tape for subsequent 
processing by the EGSA program itself. The latter is composed of a main 
program and several subprograms whose linkage is shown in the macro-flow 
diagram in Figure 3. In the discussion below we identify the purpose of 
each of these programs while detailed information is presented in the 
appendix. This was included in this report for those readers who may be 
interested in evaluating the detail content of other NASA imagery. 

The first major step in the EGSA main program is the input of user 
data on punched cards in format 2F10.0, 7A4, FIO.O. The following variables 
or arrays are specified for each data set to be processed: 

FINV - spatial frequency interval desired for output OTF data in 
cycles/mra 

FMAX - maximum spatial frequency of output OTF data in cycles/mm 

OPTION (7) - alpha-numeric array identifying the types of output 
data desired by user 

FCO - cutoff spatial frequency in cycles/mm used in smoothing the 
line spread 
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Figure 3 MACRO FLOW DfAGRAM FOR EGSA MAIN PROGRAM 
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The edge data is read from an unformatted input data tape and 
differentiated by Subroutine CONVLV. A special Fourier transform algorithm 
is employed in Subroutine FOURTH to compute the modulus (MTF) and phase 
of the Optical Transfer Function (OTF) . If requested by the user as an 
option a smoothed Line Spread is computed using truncation of the OTF at 
frequency "FCO." The remainder of the main program provides plotted, printed 
or punched card output data depending upon the options identified by the 
user. The options available include: plots of the individual edge data, 

line spreads, smoothed line spreads and transfer functions (MTF only); 
plots of the average line spreads and MTF; and punched values of the 
individual and average MTFs. 

Provision is made for processing multiple edge data sets under 
multiple options. The input data tape is read until a null record is 
encountered indicating the end of the current data set. The output require- 
ments are then completed for that edge data set; The input data tape is 
read a second time for the next edge data set. If one exists it is processed 
under the same user specified data and options as the previous set. If a 
second null record is encountered, however, the program recycles to read new 
user input data. If no input is provided execution is terminated. For 
a single edge data set the null records may be deleted and execution will be 
terminated by the same sequence of operations through recognition of the end 
of identified files on the data tape and the main stream input. 

In order to exercise the EGSA main program sixteen subprograms are 
required in addition to the standard FORTRAN I/O subprogram. Four of these 
programs are user supplied and their listings are included in the Appendix. 
CONVLV performs the convolution between two input functions or an input 
function and an internally specified function. These include an 11-point 
smoothing function and an 11-point differentiating function. DXSCAL computes 
the increment or scale factor (inches /point) for the x-axes of the individual 
function plots (^dges, line spreads and smoothed line spreads). FORINV 
performs an inverse Fourier transform based upon a special algorithms for 
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transforming a complex fvinction, inputed as modulus and phase, to a real out- 
put functions (i.e.OTF to Line Spread). FOURTR performs a Fourier transform 
of a real inputed fmiction providing a complex output function as modulus and 
phase (i.e. Line Spread to OTF) . The restrictions on the type of functions 
permits the use time saving algorithms in both FORINV and FOURTR without 
restricting the number of points in the input function as the more popular 
Fast Fourier Transform (FFT) algorithms does. 

Ten of the required subprograms are part of the Calspan Corporation 
Program Library; all but one are used to generate graphic output. AXIS draws 
an axis with tick marks at one inch intervals, annotates the value of the. 
variable at each of these marks, and provides an identification label. CLEAR 
fills each byte of a specified block of storage with zeros. This operation 
could be accomplished, less efficiently, by a sequence of do loops in the EGSA 
program. EFPLOT is used to close the plot output file after all plot output is 
completed. LINE actually causes a plot to be drawn, consisting of line 
segments connecting specified data point coordinates expressed in inches through 
two data arrays. MGRID provides a grid of variable dimensions and interval size 
in the X and Y directions independently. NUMBER converts a floating point 
number and draws the resulting alphanumeric characters on a plot. In the EGSA 
program it is used to place identification data and scale factors on the 
output plots. PLOT is a basic plotting routine used to control pen position. 

In the EGSA program it is used to reset the plotter origin. PLOTER is used 
to specify the plotting instnoment, lined or unlined paper and the type and 
color of pen desired. SCALE scan the elements in a data array to find maximum 
and minimum values, computes a optimum scale value based upon plot dimensions 
(inches) and then changes the array values to corresponding "plotter inches." 
SYMBOL draws selected symbols or special special characters on a plot. In the 
EGSA program it is used to title some of the plots. 

The remaining two subprograms are IBM supplied FORTRAN IV LIBRARY 
Subprograms. AMAXl locates the largest value in an array of real variables. 

SQRT is a mathematical function that computes the square root of a real 
argument . 
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Listings of the first group of subprograms are included in the 
appendix. The subprograms in the remaining two groups are not included since 
the user is likely to have access to the same or equivalent programs within 
his computing facility. The appendix also includes examples of the user 
supplied input data cards, description of the record sequence for the 
unformatted input data tape and examples of the output data generated. 

2.3 Tone Quality Evaluation Methods - Tone quality is related to 
density differences of large objects in an image. As the object area decreases 
in size the OTF reduces the density differences achieved between the object and 
background. Tone quality is expressed by the net D-Log E response curve which 
determines the density difference between a large object and the background 
as illustrated by the sketch in Figure 4. If we assume that apparent scene 



Figure 4 EFFECT OF D-LOG E RESPONSE CURVE UPON OBJECT CONTRAST 
AND TONE QUALITY 
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brightness is such that the difference in log exposure between the object and 
background is 

^3 ( 1 ) 

then the average exposure level (= ) influences the density difference. 

2 - 

This average exposure level includes the camera parameters (f-number, shutter 
speed, lens transmission, etc.) as well as the exposure settings used to copy 
the original film. If the level is too high or too low the density difference 
(i.e. cfi or indicated in the Figure) is less than achieved by 

exploiting the middle portion of the D-Log E response curve. This portion is 
nearly linear and is customarily represented by a straight line equation, namely 

( 2 ) 

where is the net film gamma and represents a constant gain between log 
exposure and density. At the higher and lower exposure levels the response 
curve departs from linearity; the gain is reduced and no longer is constant. 



and the resulting density difference, (f , in the constant gain portion of the 
curve determined by using Eqs. (2) and (4); 

(5) 


The ability of the eye to detect an object is approximately proportional to 
its density difference. In this case the density difference depends only on 
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the object -contrast and the net film gain. In the non-constant gain regions 
the density difference is reduced and also depends upon the exposure levels. 
Therefore ah appropriate measure of tone quality is the gain achieved as a 
function of object contrast and average exposure level, namely 

This quality parameter was calculated in this study using a spline interpolation 
between the measured step wedge densities (sensitometric calibration). The 
reader should be made aware that the step wedge densities and resulting response 
curve can depend upon the type of instrument used in measurement. In general 
macro-densities measured by a dens itonieter: (such as a MacBeth Model TD102) are 
not always equal to micrb-densities measured by a microdensitometer. The first 
instrument measures "diffuse" density while the latter generally measures 
"specular" density. The data in this study were obtained using a microden- 
sitometer and our tone quality evaluation is more appropriate to the use of 
magnifications (i. e. microscope) when viewing the image rather than unaided 
viewing. 

Calculations of gain factors, Eq. (6), were made at three levels of 
object/background contrast indicated in Table 2. The high contrast object 
would represent shadow/sunlight boundaries in Apollo imagery; the mid-contrast, 
albedo variations associated with bright rays; while the low contrast represents 
albedo changes caused by 1% reflectance variation about a 10% background 
reflectance. 


Table 2 

OBJECT CONTRASTS USED IN TONE QUALITY EVALUATION 


CONTRAST 

LEVEL 

VEr 

C 

D 

TYPICAL LUNAR OBJECT 

HIGH 

6:1 

5.0 

0.778 

SUNLIGHT/SHADOW BOUNDARY 

MID 

1.3:1 

0.30 

0.114 

BRIGHT RAY 

LOW 

1.1:1 

0.10 

0.041 

1% ALBEDO CHANGE 
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3. EVALUATION OF APOLLO 15 PANORAMIC CAMERA V/H ANOMALY 

The Panoramic Camera flown on several Apollo missions was a modified 
Itek KA-80A camera. During operation its lens system is rotated about the 
velocity axis of the Command Module at a rate related to the apparent ground 
speed. When a photograph was taken a capping shutter opened while the lens 
scanned a 108° arc below the spacecraft producing a frame about 45 inches 
long. A variable width slit located next to the film acted as a scanning 
shutter. Its velocity (scanning rate) and width determined the instantaneous 
exposure time. Image motion compensation (IMC) is provided by a complex 
electro-mechanical subsystem. Forward motion is compensated by a gimbal 
that rocks a lens group to stabilize the image. A V/H (velocity/height) sensor 
measures the apparent ground motion to determine the rate of rocking needed. 

This sensor also controls the speed of lens system rotation and consequently 
the rate of film advance which is perpendicular to the flight direction. A 
second sensor, an exposure meter, measures the brightness of the lunar 
scene which is combined with the V/H data to control the width of the 
scanning slit. 

During the Apollo 15 mission intermittent operation of the V/H sensor 
produced an anomaly in a number of the panoramic photographs. One of the tasks 
undertaken was to assess the effect of these anomalies on the detail content. 

In the photographs examined they re$ulted in (1) a loss of resolution due to image 
blur from uncompensated motion and (2) exposure differences within the frame. 

The latter effect is clearly visible in Figure 5 which shows Frames 316 and 
360 evaluated in this task. These frames were taken at low sun elevation angle 
and selected to provide shadow-to- sunlight edges for evaluation. The exposure 
differences are thought to be a direct consequence of the low sun angle. 

Telemetry data shows that the camera exposure meter estimated the average 
lunar surface brightness to be about 200 foot-Lamberts for Frame 316 and 
less than 100 foot-Lamberts for Frame 360. In both cases the telemetry data 
reproduced in Figure 6(c) shows that the focal plane or scanning slit remained 
at maximum width (0.3 inches) during the complete exposure of both frames 
due to the low surface brightness. Thus the abrupt change in the V/H rate 


13 




14 


Figure 5 TYPICAL APPOLLO 15 PANORAMIC PHOTOGRAPHS WITH V/H ANOMALIES 



MRAD/SEC MRAD/SEC 


FRAME 316 


FRAME 360 



(a) 


CAPPING SHUTTER 
PULSE AND 
FMC RATE 


(b) 

V/H COMMAND 
RATE 


(c) 

SLIT WIDTH 


_r 



Figure 6 STRIP CHART RECORDS OF TELEMETRY DATA 
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indicated in 6(b) produced a corresponding change in lens rotation and the 
velocity of the slit relative to the film without a compensating change in 
slit width. This resulted in the exposure differences. In frames taken at 
higher sun angles or brighter lunar surfaces the slit width could be adjusted 
permitting reduction or elimination of these exposure changes. 

The abrupt change in the commanded V/H rate was produced by the V/H 
sensor signal "going into" or "coming out of" saturation. When the signal 
exceeded 19. S mrad/sec (established as the maximum rate that could be 
encountered in the Apollo missions) the commanded V/H rate was automatically 
reset to 13.0 mrad/sec (a nominal value) producing the abrupt changes shown 
in Figure 6(b). When these changes occurred while the capping shutter 
was open. Figure 6(a), the detail content of the photography was reduced. 

To estimate the amount of residual motion blur produced by this 
V/H rate anomaly, microdensitometer traces were made of the shadow-to- 
sunlight edges on both sides of the abrupt change in commanded V/H rate. 

From the telemetry data presented in Figure 6 we see that the V/H sensor 
signal reached the saturation value and the V/H command was set to the 
nominal value for both of the frames evaluated. The first portion of 
Frame 316 is dark and has a commanded V/H rate varying from 17.6 to 19.3 

\ 

mrad/sec and the latter portion is light with a V/H rate of 13.0 mrad/sec. 
Similar conditions exist for Frame 360 except that the commanded V/H rate 
is nearly constant equal to 19.2 mrad/sec in the dark portion of the frame. 

None of the commanded V/H rates resulted in full compensation of the image 
motion. The residual blur was measured using the EGSA prograim and also 
calculated directly from telemetry and auxiliary data. The results of both 
of these estimates are compared in this section. 

3.1 Residual Blur Estimation - Telemetry and auxiliary data was used 
to estimate the residual blur in both the dark and light portions of each of 
the frames evaluated. This was achieved by calculating the difference 
between the uncompensated angular blur and the amount of compensation provided. 
In performing these calculatibns , estimates of the angular blur rates and 
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exposure times were made using the available data sources and measurements 
of the photographs themselves. We discuss the calculations of the uncompensated 
angular blur rate, the compensation rate, the exposure time and the residual 
blur separately in the subsections below. 


3.1.1 Uncompensated Angular Blur Rate - In order to compute the 
angular blur rate due to the velocity of the Command Module provision must 
be made for the curvature of the lunar surface. In a previous study 
expressions for the angular blur near the principal point in vertical and 
forward or aft oblique photography were developed. The expression for the 
ground blur at the lunar surface is 




Vte 






(7) 


where V is the velocity of the Command Module, is the exposure time, 
h is the altitude of the Command Module and is the nominal lunar 

radius (= 1738.1 kilometers). As the altitude decreases, this expression 
is asymptotic to the more familiar expression in which the blur is proportional 
to the product of the platform velocity and the exposure time. For Apollo 
orbital photography, the reduction in the blur due to the interaction of the 
Command Module altitude and the lunar surface curvature is approximately 6%. 


In our earlier study, all assessments were made near the principal 
point in each photograph and the theoretical expressions developed are only 
valid under these conditions. In order to evaluate the Apollo 15 V/H nominally, 
our assessments had to be made away from the principal point in order to 
evaluate the MTF and residual blur on both sides of the V/H anomaly. Figure 7 
shows a sketch where a point in a forward oblique pan photograph is under 
evaluation. The point is specified by a pitch angle, p , and a roll angle, zi/" . 
In general the pitch angle would be the sum of the camera pitch angle plus 
an additional term due to the location of the point across the width of the 
photograph. In the present case however, all assessments were made near the 
optical axis ground track where the additional term is negligible. The pitch 




Figure 7 PAN PHOTO GEOMETRY 


angle will result in a decrease in the apparent ground blur due to fore- 
shortening as discussed below. The roll angle, on the other hand, being 
perpendicular to the Command Module ground track does not produce any fore- 
shortening effects, however, the ground blur distance will be decreased from 
that expressed by Eq. (7) due to the fact that the point under assessment is 
away from the principal point by an amount expressed through the roll angle. 
For a vertical photograph an additional reduction factor, cose , must be 
included in Eq. (7) so' that 


where 



\/t^ cos O 


1 ^ hi 


m 


(8a) 


Q= stn'^ 


stn Ticr 
_ 


(8b) 


and is the slant range from the Command Module to the point under 
evaluation. The altitude of the Command Module is sufficiently small 
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compared to the lunar radius in this case so that this reduction factor is 
negligible in vertical panoramic photography. For oblique photography, the 
pitch and roll angles interact producing a more complicated expression for 

r 

the reduction factor, however, its effect remains negligible. 


The next step is the computation of the apparent ground blur by 
including the effects of foreshortening on the actual ground blur given by 
Eq. (7). Both the roll and the pitch angle can contribute to the foreshortening 
effect through the tilt angle given by 


Cos t= cos <f) cos ur 


C9) 


When the pitch angle is zero, however, the roll angle produces no foreshortening 
effect since all of the ground blur is perpendicular to the line of sight. 

The foreshortening factor, cos0 , was included in the expressions presented 
in our previous study. In the presence of a roll angle the foreshortening 
factor is cos ^ as presented by Eq. (9). However, this factor should only be 
applied to the portion of the blur in the viewing direction. This requires 
knowledge of the azimuth angle between the direction of flight and the 
direction of tilt. Kawachi presents the required transformation equations 
permitting us to write the apparent blur as 


\/t^ fsin^oL /• cos^cL cos^t) 
m 

where the azimuth angle, oc , is given by. 


(10a) 


ce:= tdzi ^ 


sin ur 

t<3.n 0 


(lOb) 


For the values of roll angles which we are considering, however, this 
expression can be approximated by the original on-axis equation, namely 

\/t^ cos (p 

A “ -jTj;; (10=) 

m 
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This expression is a valid approximation because of two opposing effects; 
although the foreshortening factor is smaller, i.e. cos i compared to Cos0 , 
it is only applied to a smaller portion of the actual ground blur. These 
two facts compensate for each other, allowing the foreshortening effect to 
be approximated by cos(p with an error less than 0.1%. 


The foreshortening effects, including the geometrical equations 
developed by Kawachi, assume a flat planetary surface. As a result of this 
assumption, the apparent ground blur will be slightly overestimated. Figure 8a 
shows the percentage error by which the blur is overestimated for two 
representative Command Module altitudes. For the present case, the error is 
less than 0.5% since the pitch angle will be either 0° or 12.5° determined 
by the. mode of camera operation. 

The next step in the calculation of the uncompensated angular blur 
rate is the conversion of the apparent ground blur to an apparent angular 
blur by the use of the slant range from the Command Module to the ground . 

point under evaluation. In our earlier study, the expression for slant 
range employed was 



/? 

cost 


Cll) 


which is based upon the assumption of a flat planetary surface. The actual 
slant range is larger due to the curvature of the lunar surface and is given 

Figure 8b shows the error produced in the estimation of slant range by assuming 
a flat lunar surface. The net effect of the assumption is to overestimate 
the angular blur similar to the effect on foreshortening. In this case, 
however, errors on the order of 1% are expected since tilt angles can become 
as^fhigh as 30°. Finally, the uncompensated angular blur rate is determined 
by taking the apparent ground blur of Eq. (10c) and dividing by the slant 
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(a) FORESHORTENING EFFECT ERROR 
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(b) SLANT RANGE ERROR 

Figure 8 ERRORS RESULTING FROM THE ASSUMPTION OF A FLAT LUNAR SURFACE 
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range expressed by Eq. (12) and the exposure time, • The uncompensated 

angular blur rate, yS , is 

• // cos (p 

( 1 fh/rj ff/) cos t f ^sin J (13) 

This expression was employed with the numerical values shown in Table 3 to 
calculate the blur rate for both the dark and light portions of the frames 
under evaluation. The values of the Command Module velocity and altitude 
and the pitch angle were obtained from auxiliary data furnished by Mr. K. 
Hancock at the NASA Johnson Space Center. The roll angles were estimated by 
measuring the distance of the dark and light evaluation regions from the 
center of the frame which is directly proportional to the roll angle. The 
tilt angles are subsequently calculated from Eq. (9). 


Table 3 

NUMERICAL VALUES USED TO COMPUTE 
UNCOMPENSATED ANGULAR BLUR RATE 


VARIABLE 

SYMBOL 

VALUES 

FRAME #316 

FRAME #360 

CM VELOCITY 

V 

1.6279 km/sec 

1.6305 km/sec 

CM ALTITUDE 

h 

109.22 km 

106.14 km 

LUNAR RADIUS 

m 

1738.1 km 

1738.1 km 

PITCH ANGLE 

0 

12.45° 

0 

ROLL ANGLES 

Wd 

20.4° 

6.0° 

Wl 

25.2° 

10.8° 

TILT ANGLES 

^D 

23.8° 

6.0° 

*L 

27.9° 

10.8° 


The following values resulted; Frame 316 

mrdd/sec (14a) 

/d^ = n-9? mrdi/sec 
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and for Frame 360 


P^= 14 . 3 9 mvdd/ Sec 
/3^ = 14. 20 Tnraj/sec 


(14b) 


An expression equivalent to Eq. (13) could also have been developed 

(3 4) 

by employing the work of Kawachi^ ' . He derives expressions for the 

f4) 

image velocity for an oblique frame camera and notes ^ ^ that "the image 
velocity for the panoramic photograph is the same as for a frame photograph 
tangent to it at the point in question." If we are interested in evaluating 
the velocity near the center of the width of the panoramic frame we can set 
the transformed film coordinates equal to zero, , in Eqs. (7) 

and (8^ of Ref. 4 and obtain the net image velocity as 




/// 


COS t \_sm ^cc f cos cos 


^ 4 . 


(15) 


where /4 is the focal length of the camera and od is defined by Eq. (10b) . 
The uncompensated angular blur rate is determined by dividing by the product 
of the focal length and the factor included previously in Eq. (7). 

This factor is the principal correction for the curvature of the lunar 
surface. The resulting expression can also be obtained from Eqs ..(10a) and 
(LI) and is therefore equivalent to Eq. (13) except for the neglection of the 
effect of surface curvature upon slant range. 


3.1.2 Blur Compensation Rate -The next step in computing the 

residual blur is the determination of the compensation rate applied by the 

( 2 ) 

camera. In the case of a panoramic camera the blur compensation rate is 
proportional to the instantaneous rotation angle of the lens, namely. 


/5c 



(16) 


This reflects the decrease in angular blur due to the slant range assuming 
a flat planetary surface. It is not known whether. the panoramic camera IMG 
computer was reprogrammed for the more complex effect of the slant range due 


23 



to the cin-vature of the lunar surface (c.f. Eq. (12)). If so, the compensation 
rate would be 




/ 


H/r 


^2^7'^/^J CoSli^- 3tn^W ^ (17) 


where can be assumed to be constant. Referring back to Figure 8b, 

we see that the effect is less than 1% for the rotation or roll angles we 
are considering. However, near the edge of the frame where the rotation angle 
approaches 54° the error would be as large as 7%. Conversely, the uncompensated 
blur rate is considerably less in this portion of the panoramic frame and the 
need for an accurate compensation rate diminished. Consequently, it is assumed 
that such a modification of the IMC computer was not undertaken since it would 
offer little advantage. Using Eq. (16), the roll angles presented in Table 3 
and the V/h command values from the telemetry data presented in Figure 6, 
the compensation estimates were estimated for Frame 316 to be 


and for Frame 360 


- 18. 00 /sec 

= 11-76 mrad/sec 


(18a) 


/6^^=iaoS ?nZHd/3ec 

/d„ ^ /2-77 

We note that these results are only an approximation since the 
compensation rate is proportional to the instantaneous rotation angle and 
a given Point in the panoramic frame will receive a varying amount of image 
motion compensation due to the finite width of the scanning slit. Consequently, 
Eq. (16) should be integrated over the exposure time of the point under 
evaluation to determine an average compensation rate, namely 




( 19 ) 
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This procedure was executed and the resulting average compensation rates 
were less than those presented in Eqs. (18a and 18b) by approximately 0.6%. 

In this case, therefore, the dynamic effect of the motion compensation can be 
neglected. 


3.1.3 Exposure Time Estimation - Three different methods were used 
to estimate the exposure time of the dark or light portions of each frame 
under evaluation. As noted previously, the abrupt change in the average 
density across the V/H anomaly was attributed to a change in exposure time 
between these two portions of the frames. One method for determining the 
exposure times is based upon this fact, that is, the density difference is 
related to the ratio of the exposure times by the expression 

As. = (20) 

where cT is the density difference and 7" is the gain of the D-Log E response 
curve. The length of time that the capping shutter is open, is equal to 

the sum of the times required to , expose each portion of the frame, that is, 






A 

/I 




!/ / 


( 21 ) 


where is the length of the dark portion of the exposed frame, is the 
length of the light portion of the exposed frame, and 1^ , are the 
velocities of the scanning slit in the respective portions of the frame. 
Since the width of the scanning slit, > is constant these velocities 
are inversely proportional to the exposure time and consequently we can 
equate the variable of Eq. (20) to the ratio of the velocities, namely 



( 22 ) 


Combining Eqs. (21) and (22) and eliminating the velocity of the scanning slit 
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in favor of the ratio of the slit width to the exposure time, we can derive 
an equation for the exposure in the light portion of the frame 







4 


as well as an expression for exposure in the dark portion 


(23a) 




(23b) 


These equations were employed with the numerical values shown in Table 4 to 
calculate the exposure times for both the dark and light portions of the frames 
under evaluation. The density difference was measured directly from the 
frames under evaluation while the gain of the response curve was obtained 
from Figure 31 of Reference 1 using the curve labeled "4th GEN POS (4P)", 

The capping shutter time and the scanning slit width were obtained from 
telemetry data shown previously in Figure 6. The values of the exposure 
time -Obtained using these expression are presented in Table S along with 
error estimates derived from first order analyses. 


Table 4 

NUMERICAL VALUES USED TO ESTIMATE EXPOSURE TIMES 


VARIABLE (UNITS) 

SYMBOL 

VALUE 

FRAME # 316 

FRAME it 360 

DARK-LIGHT DENSITY DIFFERENCE 

5 

0.42 

0.36 

RESPONSE CURVE GAIN 

7 

-2.0 

-2.0 

EXPOSURE TIME RATIO 


0.61660 

0.66069 

SECTION LENGTHS 
(INCHES) 

•-D 

32 

26 

k 

13 

19 

CAPPING SHUTTER TIME (SECONDS) 


1.5 

3.1 

SCANNING SLIT WIDTH (INCHES) 

Ws 

0.3 

0.3 
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Tjdile 5 

EXPOSURE TIME ESTIMATES 


SOURCE 

FRAME ;l'316 

FRAME #360 

DARK 

LIGHT 

DARK 

LIGHT 

EOS (23a) AND (23b) 

8.5 msec 
± 9.3% 

13.7 msec 
+ 7.7% 

17.0 msec 
± 5.3% 

25.7 msec 
± 4.2% 

MOLEBURG EQ. 
(25) 

9.2 msec 
± 2.9% 

13.6 msec 
± 2.6% 

17.6 msec 
+ 2.3% 

26.0 msec 
+ 2.6% 

ADJACENT FRAME 

N/A 

14.0 msec 
± 5 . 0 % 

N/A 

25.3 msec 
± 3.2% 


The second method used to estimate the exposure times is based upon 
an approximate formula which was furnished by Mr. B. Moleburg at the NASA 
Johnson Space Center, namely 

where ^ =0.59 for the stereo mode and 1.127 for the mono mode. The V/H 
factor in the denominator of this expression is the command V/H rate obtained 
from telemetry data; 19.2 milliradians/second for the dark portion and 13.0 
milliradians/second for the light portion. The exposure times obtained using 
this expression are also presented in Table 5 with their associated error 
estimates . 

The final method used to estimate the exposure time is only applicable 
to the light portion of each frame under evaluation. In both cases, the 
adjacent frame had an exposure level equivalent to the light portion of the 
frames under evaluation and consequently the exposure time in the light 
portion is equal to the exposure time of either adjacent frame. The equation 
for calculating the .exposure time for a frame with a constant V/H command rate 
is obtained by setting!^ =1, =0 and =45 inches in Eq. (23a), 

namely 
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( 25 ) 



The results obtained using this -method are also included in Table 5, and we 
note that the results from all three methods are in excellent agreement. 

In subsequent evaluations, results obtained from Eqs. (23a) and (23b) are 
employed. 


3.1.4 Residual Blur - It is a straightforward procedure to combine 
the results derived above in. order to determine the residual motion blur at 
the film plane. The equation for the residual image blur is 


V|/^- A| 4^ 


(26) 


where f is the focal length of the camera (=610 millimeters). This blur 
produces a Modulation Transfer Function described by the equation 




Sin f7rv/Po) 


(27) 


which degrades the nominal fine detail content of the photograph. The 
parameter t?o is simply the reciprocal of the residual image blur in cycles/mm 
Values obtained for both the residual image blur and the "cutoff" spatial 
frequency are presented in Table 6. 


Table 6 

RESIDUAL IMAGE BLUR 



FRAME 

FRAME #360 

VARIABLE (UNITS) 

DARK 

LIGHT 

DARK 

LIGHT 

RESIDUAL IMAGE BLUR 

.0288 

.0020 

.0224 

.0488 

(mm) 

+ 12% 

±44% 

±7% 

± 17% 

CUTOFF FREQUENCY 
(CYCLES/mm) 

34.7 ± 4.2 

512 ±230 

44.5 + 3.1 

20.5 + 3.5 
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The expected value of the net WTF can be determined by multiplying 
the values obtained from Eq. (27) by the product of the MTF for the flight 
and subsequent copy films and the camera lens. Unfortunately, an MTF for 
the panoramic camera lens was not available, however, such lenses are noted 
for their high quality performance and it is very likely that the net MTF is 
dominated by the performance of the flight and copy films and any residual 
blur. Consequently, in determining the expected MTF the lens MTF was 
assumed to be 1.0 at all spatial frequencies evaluated. 

A final comment is required concerning the precision of the 
approximations used in the expressions for the uncompensated angular blur 
rate and the compensation rate. In both cases only those approximations 
whose errors were less than 1% were retained. However, since we are employing 
difference between these two quantities in Eq. (26) a small error can be 
magnified when a small difference occurs. Eq. (13) for the uncompensated 
angular blur rate employs two approximations for the inclusion of foreshortening 
effects; the assumption of a flat planetary surface and equivalence between 
the off-axis and on-axis reduction factors. For vertical photography, these 
approximations introduce no errors. However, for oblique photography, the 
net effect is to over estimate the blur rate by about 0.5%. Eq. (16) for 
the blur compensation rate neglects the dynamic effects of motion compensation 
and over estimates by 0.6% compared to the average compensation rate. This 
occurs for both the oblique and vertical photography. In the case of 
oblique photography, since both the blur rate and compensation rate are 
over estimated by about the same percentage, the effect upon the residual 
image blur is negligible (i.e., about 0.6%) even when there are small 
differences between the rates. In the case of the vertical photography, 
only the compensation rate is in error and this is subsequently magnified 
by the differencing operation resulting in an under estimate of the residual 
•image blur by 5% in the light portion of Frame 360 and an over estimate of 
the residual image blur by 3% in the dark portion of Frame 360. These errors 
are still less than the net accuracy of the residual blur estimates shown in 
Table 6. 
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3.2 Modulation Transfer Fiinction Estimation - As noted previously 
in Section 2, a net MTF was estimated for each photo under evaluation. This 
was accomplished by tracing shadow-to-sunlit edges interior to craters using 
a microdensitometer and subsequently processing this edge data using the 
Edge Gradient Spectral Analysis (EGSA) software package described in Section 2.2. 

A Mann microdensitometer was used to obtain the required edge data. 

For this purpose the microdensitometer was used in an improved but unconventional 
mode of operation. A 41//m by 41 /im source aperture was employed with a 0.4 N.A. 
influx objective. This offered a magnification of approximately 40x producing 
a 1.0 ><m square geometric image at the film plane. This image is subsequently 
broadened by the diffraction effects of the objective whose airy disc has a 
diameter of 1.5// m. Consequently, the image of the source aperture will be a 
gaussian-like light distribution with a diameter of several micrometers. A 
0.25 N.A. objective was used for the efflux optics and the subsequent sensor 
aperture chosen to. minimize stray light and not further define the scanning 
aperture. In this mode of operation the distorted image of the source 
aperture is the scanning aperture and the spatial frequency response of the 
raicrodensitoraeter determined accordingly. The over filling of the efflux 
optics by the influx optics offers an improvement in instrument performance. 

Edge data was collected in both the dark and light portions of each 
of the frames about 1 inch on either side of the suspected V/H discontinuity. 

This edge data was then hand smoothed and digitized for processing by the 
EGSA software. The average MTFs resulting from these measurements are 
presented in Figures 9 and 10 along with the expected MTF determined by the 
procedures described in Section 3.1.4. In estimating the average MTF through 
modulus averaging, 10 edges were measured in the dark portion of Frame 316 
while 7 edges were measured in the light portion. Similarly, 11 edges were 
measured in the dark portion of Frame 360 and 10 edges in the light portion. 


■o . 
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MODULATION TRANSFER FUNCTIONS 



SPATIAL FREQUENCY (cycles/mm) 


Figure 9 MEASURED AND EXPECTED MTF's FOR APOLLO 15 FRAME 316 WITH 
V/H ANOMOLY 
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1.0 



Figure 10 MEASURED AND EXPECTED MTF's FOR APOLLO 15 FRAME 360 WITH 
V/H ANOMOLY 
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Also shown in the Figures is the net film MTF limit for the 4th Generation 
Positive copies employed in the evaluation. We note that there is excellent 
agreement between the expected and measured MTF values in all four instances. 
For the light portion of Frame 316 the measured and expected MTF values are 
essentially identical and equal to the net film MTF limit over the range of 
spatial frequencies shown. This is a result of the nearly perfect compensation 
provided in this case. In all the other cases the compensation was poor and 
the resolution reduced by 1/3 (to about 6 meters on the lunar surface) . 

Other low sun angle photography known to have V/H discontinuities 
include Frames 8852 (Rev. 4) and 344 (Rev. 72). When such a discontinuity 
did not occur the V/H command signal was either at the nominal value (13.0 
mrad/sec) similar to the light portions of the frames evaluated or varying 
about an anomalously high value similar to the dark portions. Telemetry 
data indicates that the latter situation occurred at least for Frames 345, 

346 and 347 (Rev. 72). Only when it is equal to the nominal value can the 
full "resolution" potential of the panoramic camera be realized; but this 
may not always occur (refer to the MTF for the light portion of Frame 360) . 

The user of Apollo 15 panoramic data should refer to the telemetry data 
to determine the V/H command situation for the frames he is exploiting and 
account for the effects of poorer resolution if necessary. 



4. COMPARISON AMONG MULTIPLE GENERATIONS OF APOLLO 16 AND 17 PHOTOGRAPHY 

One of the specific tasks accomplished during this study was an 
assessment of the detailed content of multiple generations of Apollo 16 and 
17 photography. The purpose of this assessment was to determine the 
degradation, if any, in the quality of the photographic products furnished 
to variQus users. The original flight film is copied to produce a Master 
Positive (2P) of Direct Negative (2N) copy.; These are the earliest generations 
which a user may exploit and only a limited number of such copies are generated 
The Master Positive is subsequently copied^ twice to produce a 4th Generation 
Positive (4 P) which most investigators are likely to use. The Direct Negative 
is also copied once to produce a Third Generation Positive (3P) which is also 
available to investigators. This latter sequence eliminates one copy step 
in the hope of improving product quality or reducing production time. Unfor- 
tunately the characteristics of the reversal film required to produce the 
2N copy may limit the detail content and no actual improvement in product 
quality obtained. The purpose of this task was. to assess the effect of both 
reproduction sequences on the detail content of Apollo 16 and 17 photography. 

An inter-comparison of the total quality and the fine detail content 
was made for Apollo 16 panoramic photography, Apollo 17 panoramic photo- 
graphy and Apollo 17 metric photography. The frames selected and subsequently 
evaluated were shown previously in Table 1. Unfortunately 3P and 4P copies 
of the same type of photography were not available at the Lunar Science 
Institute thus preventing a direct comparison of the most likely user products. 
Instead a comparison of proceeding generations (2N and 2P) was made. The 
results of the evaluations are discussed in detail in the subsections below. 

4.1 Product Tonal Quality - The methodology used to assess the tonal 
quality of the Apollo 16 and 17 photography was discussed previously in 
Section 2. The most common measure of tonal quality is the Hurter-Driffield 
or D-Log E response curve. It specifies how the apparent scene brightness 
is reproduced as density in the photograph. The D-Log E response curve was 
measured for each frame evaluated. The results obtained are presented in 
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Figures 11 through 13. These curves were generated from raicrodensitometer 
scans of step wedges placed on the original flight film and subsequently 
copied during each duplication step. Consequently, the relative exposure 
refers to the apparent brightness of scene elements at the original flight 
film. The exposure of the first step was arbitrarily assigned a value of 
1.0. Since there is an increment of 0.15 in the log exposure between the 
21 steps, the last step has a relative exposure of 1000. 

A software program was written that calculates the gain factors 
expressed by Eq. (6) using step wedge densities (sensitometric calibration) 
supplied as input. The input values were obtained from the curves shown in 
Figures 10 through 13 created by visually passing a smooth curve through the 
actual measured step wedge densities. Spline interpolation functions are 
used to compute density values for exposure levels between the individual 
steps. These interpolation functions yield results equivalent to forcing 
a flexible straight line (spline) through the data points and tend to avoid 
the "waviness” often encountered by other curve fitting techniques. 

The tone quality evaluation program uses the resulting coefficients 
of the spline regression to estimate the density difference between the 
image and the background for various contrast levels identified previously 
in Table 2. The average exposure levels, i.e. the average of the object 
and background exposure values, were set equal to the relative exposure 
levels of each of the steps in the step tablet. The density difference 
at each average exposure level is divided by A (c.f. Table 2) to determine 
the gain. 

The gain as a function of the average exposure level was identical 
for the low and mid-contrast objects for all the cases examined. The gain 
functions for the second generation photography are presented in Figures 14 
through 16. In describing tonal quality the term "dynamic range" is frequently 
employed. This is loosely defined as the exposure region corresponding to the 
linear portion of the D-Log E curve. In this study we defined the dynamic 
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Figure 11 SENSITOMETRIC CALIBRATION - APOLLO 16 PANORAMIC PHOTOGRAPHY 
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MICRODENSITY MICRODENSITY 



4 6 8 10 20 40 60 100 200 400 600 1000 

RELATIVE EXPOSURE 
(a) FRAME 1633 



4 6 8 10 20 40 60 100 200 400 600 1000 

RELATIVE EXPOSURE 
lb) FRAME 3111 


Figure 12 SENSITOMETRIC CALIBRATION - APOLLO 17 PANORAMIC PHOTOGRAPHY 
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MICRODENSITY MICRODENSITY 



RELATIVE EXPOSURE 
(b) FRAME 2923 


Figure 13 SENSITOMETRIC CALIBRATION - APOLLO 17 METRIC PHOTOGRAPHY 
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PANORAMIC PHOTOGRAPHY 





2P COPY 
2N COPY 
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PANORAMIC PHOTOGRAPHY 
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METRIC PHOTOGRAPHY 





range more quantitatively as the exposure interval between the points where 
the gain functions dropped to 50% of their peak value. The average value for 
the peak gain and the dynamic range for each type of photograph evaluated 
is presented in Table 7. It is clear from examination of Figures 14 and 15 
that both the 2P and 2N photography have about the' same peak gain but the 2P 
photography has a much greater dynamic range. This is also evident from the 
numerical values presented in Table 7. We can conclude, therefore, that the tone 
quality of the Master Positive is superior to the Direct Negative for the 
panoramic photography. In the case of the metric photography, it is not clear 
from Figure 16 that a difference in dynamic range exists and the values in 
Table 7 show an equivalence of the dynamic range. However, the peak gain 
of the Master Positive photography is clearly greater. The roll containing 
Frame 180 was, in fact, processed to yield a "high gamma" copy. In this case 
we conclude that the Master Positive is superior due to the peak gain rather 
than any difference in dynamic range. 


Table 7 

TONE QUALITY COMPARISON 
(LOW AND MID CONTRAST OBJECTS) 


PHOTO TYPE 

— 

PEAK 

GAIN 

DYNAMIC 

RANGE 

(dB) 

APOLLO 16 

■9 

1.8 

9.9 

PAN 

KB 

1.9 

15.5 


3P 

1.8 

9.6 



2.1 

7.5 

APOLLO 17 

Bl 

2.2 

11.1 

PAN 


3P 

Z4 

6.5 


Bl 

1.9 

8.3 

APOLLO 17 

H 

Z7 

8.4 

ME 1 KIC 


4P 

2.1 

7.5 
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Similar measurements were made. on the 3P and 4P copies where they were 
available. These data were not included in Figures 14 through 16 for the sake 
of clarity but are presented in Table 7. The quality of these products generally 
correlates with the corresponding Second Generation photography. There is 
little change in the peak gain between these generations, but a definite 
decrease in the dynamic range for the higher generation products. An exception 
is the 2P to 4P reproduction of the Apollo 17 metric photograph where an 
apparent drop in the peak gain is noted. Consequently, we cannot infer that 
the tonal quality differences present in the Second Generations are preserved 
in the 3P and 4P user products for the metric photography. 

The gain functions for a high contrast object C^.g. sunlight-shadow 
boundary) have lower peak gains and greater dynamic range (c.f. Table 8) but 
support the conclusions reached above. 


Table 8 

TONE QUALITY COMPARISON 
(HIGH CONTRAST OBJECTS) 


PHOTO TYPE 

PEAK 

GAIN 

DYNAMIC 

RANGE 

(dB) 

APOLLO 16 


1.5 

13.0 

PAN 

n 

1.7 

> 15.7 


3P 

1.5 

12.0 

APOLLO 17 


1.7 

9.7 

PAN 

' 

n 

2.0 

12.3 


3P 

1.9 

9.2 

APOLLO 17 


1.6 

10.7 

METRIC 

n 

2.6 

10.7 


4P 

1.6 

10.7 


4.2 Fine Detail Content - MTF measurements were made in order to 
compare the detail content of the various generations of photography evaluated. 
These measurements were made by scanning shadow-to-sunlit edges near the center 
of each frame using a microdensitometer. The mode in which the instrument 
was operated was described previously in Section 3.2 in connection with 
similar measurements on Apollo 15 panoramic photography. The edge data 
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obtained were processed using the EGSA software as discussed previously in 
Section 2.2. Table 9 shows the amount of edge data collected and processed 
for each of the photographs evaluated. The MTFs measured for the Apollo 17 
metric photography are presented in Figure 17. Included in the Figure is the 
NfTF limit for the film alone. The significance of the departure of the measured 
MTFs from the film limit was assessed using the 95% confidence bands calculated 
by the EGSA main program (see the Appendix). The systematic affect at the low 
spatial frequencies (0-30 cycles/mm) where the measured functions exceed the 
limit is a result of the manual smoothing and the noise level. The top and 
bottom of the edge are masked by the noise and the observer tends to employ 
small radii of curvature thereby truncating the tails of the spread functions 
and producing this effect. Consequently the confidence bands are examined only 
above 30 cycles/mm to assess significance. Both the Direct Negative (a) and 
Master Positive (b) MTF's are not significantly lower than the film limit 
except for the 2 P copy of Frame 180. The 4th Generation Positive MTF's (c) 
are both significantly less than the film imposed limit. In these three cases 
the estimated ground resolution is about 30 meters (versus a nominal 20 meter 
resolution). Since half of the copies have nominal resolution the loss in 
fine detail must have occurred during the reproduction processes. 


Tables 

EDGE DATA COLLECTED AND PROPOSED FOR 
APOLLO 16 AND 17 PHOTOGRAPHY 



DIRECT 

NEGATIVE 

MASTER 

POSITIVE 

3rd GEN. 
POSITIVE 

4th GEN. 
POSITIVE 


8 

7 

' 


10 

2293 

8 

8 

- 

8 - 

1633 

10 

9 

6 

- 

■BDH 

6 

6 

3 

- 

IQgllll 

5 

10 

9 

- 


10 

8 

7 

- 
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Figure 17 MODULATION TRANSFER FUNCTIONS - APOLLO 17 METRIC PHOTOGRAPHY 







Figure 18 presents the MTF results for the Apollo 17 panoramic 
photography. For both frames the Direct Negative, Master Positive and 3rd 
Generation Positive MTF's are not significantly different from the respective 
film limits. This photography has the maximum possible fine detail content 
or about 2 meter "resolution." 

Figure 19 shows the MTFs obtained for the Apollo 16 panoramic 
photography evaluated. In this case all of the copies examined showed a 
significantly lower response than that imposed by the photographic films 
involved. For the Direct Negative copies the ground resolution element is 
estimated to be 4 meters while the Master Positive and 3rd Generation Positive 
copies have a 3 meter resolution. Whether this reduced fine detail content 
is present in the original flight film it is unknown since that photography 
was not available for evaluation. We note that the Direct Negative copies 
evaluated are not necessarily those used to produce the 3rd Generation Positive 
copies (which apparently have better resolution). 

Bas'ed upon the results presented above it is concluded the reproduction, 
processes may have produced a loss in fine detail. About half the frames 
evaluated showed some loss. In these cases the ground resolution element 
increased by a factor of 1.5 to 2. 
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Figure 18 MODULATION TRANSFER FUNCTIONS - APOLLO 17 PANORAMIC PHOTOGRAPHY 
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Figure 19 MODULATION TRANSFER FUNCTIONS - APOLLO 16 PANORAMIC PHOTOGRAPHY 







5. CONCLUSIONS 


The development and application of image evaluation methods for 
assessing the detail content of Apollo Orbital Photography has been demon- 
strated. Edge analyses using shadow-to-sunlight edges interior to craters 
was successfully used to evaluate the fine detail content for Apollo 15, 16 
and 17 imagery. A method for evaluating tone quality was developed using 
a gain factor as a function of object contrast and average exposure level 
that can be related to object detectability. 

The anomalous operation of the V/H sensor during the Apollo 15 mission 
introduced several degrading effects into the panoramic imagery; a loss in 
resolution due to image motion blur and exposure differences within a frame. 

The latter effect is present in those frames where an abrupt change in the 
V/H rate occurs. Two such frames were evaluated and exhibited a 0.4 density 
difference. It is estimated that this difference would be less in high sun 
angle photography although no such frames were examined. 

The fine detail content was seriously degraded by the anomalous V/H 
sensor operation. The commanded V/H rate was either anomalously high or 
equal to the pre-set nominal rate. This resulted in either over or under 
compensation and a reduction in the resolution by 1/3 (to about 6 meters on 
the lunar surface) . In one instance the nominal V/H rate was sufficiently close 
to the actual rate to permit nearly perfect, compensation. Only when the commanded 
V/H rate is equal to the .nominal value can the full resolution potential of 
the panoramic camera be realized; but this may not always occur . When the 
command rate varies about an anomalously high value a loss in resolution can be 
expected for the full frame. The user of Apollo 15 panoramic imagery should 
refer to the telemetry data to determine the V/H command situation for the 
frames he is exploiting and account for the effects of poorer resolution if 
necessary. 
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As a result of the comparison among multiple generations of Apollo 16 
and 17 photography it was concluded that the tone quality of the Master Positive 
is superior to that of the Direct Negative for the panoramic photography. This 
difference is attributed to a greater dynamic range of the 2P photography. This 
means that lunar objects have a greater contrast over a wider range of average 
exposure levels in the 2P photography compared to the 2N photography. Unfor- 
tunately 3P and 4P copies of the same type of photography were not available 
thus preventing a direct comparison of the most likely user products. However, 
the quality of these products generally correlates with the corresponding 
2nd Generation Photography, and we can infer that the tonal quality differences 
are preserved in higher generation user products. Such differences were not 
evident in the metric photography. 

MTF measurements made to compare the detail content of the various 
generations revealed that half of the imagery examined had a significant loss 
in fine detail content. These losses were found in both the metric and panoramic 
photography and are probably introduced during the reproduction process. In 
the case of the Apollo 16 panor^ic imagery all copies showed loss of fine 
detail and it was not possible to assess whether this loss was present in the 
flight imagery since it was not available for evaluation. The net effect was -an 
increase in the ground resolution element by an approximate factor of 1.5; 

3 meters for the panoramic photography and 30 meters for the metric photography. 

The success of the image evaluation methods in assessing the detail 
content of Apollo orbital photography and the importance of this detail 
content to potential users suggests that NASA should continue to develop 
the methods for applications for the future manned and unmanned spacecraft 
involved in planetary or earth exploration. 
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APPENDIX 

SOFTWARE DESCRIPTIONS 


This Appendix contains descriptions of eight programs used in the 
evaluation of the Detail Content of the Apollo Orbital Photography. The 
purpose, description of inputs, logic and processing, and a sample of the 
input and output data are included, if appropriate, for each of the programs. 

All of the programs were written using FORTRAN IV Language. Source and cross 
reference listings are included for seven of these programs. 

Seven of the programs are part of the Edge Gradient Spectral Analysis 
(EGSA) software package. This package consists of two main programs and five 
subprograms and is described in Sections 1.0 through 2.0 of this Appendix. 
Section 8.0 describes the main program used to evaluate the tone quality of 
Apollo Orbital Photography. 

1.0 Main Program DTAPE - This program is part of the Edge Gradient Spectral 

Analysis software package. It is one of two main programs in the package and . 
is used to create an input data tape for subsequent processing by the main 
EGSA program. The input to this program is assumed to be edge trace data (i.e., 
density values at uniform spatial increments) supplied on card input. Sensito- 
metric calibration (step wedge data) is also required for program execution. 
These data are used by function subprogram EXPOS to convert the density values 
to relative exposure values which are subsequently written onto an output 
data tape. No printed output data is furnished by this program. 

Fiugre 20 is a listing of program DTAPE. The function of the major 
blocks are identified by the commen^ts included in the coding. The program 
initially reads the sensitometric calibration data supplied by the user in 
card form. An example of a typical input data set is shown in Figure 21. The 
initial card of this data deck specifies the number of points, NSW, and the 
sense of the photographic image (i.e., positive or negative). The variable 
SENSE is actually a dual purpose variable. It is used to indicate the sense 
of the image; a positive value corresponds to a positive image and a negative 
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Figure 20 SOURCE AND CROSS REFERENCE LISTINGS FOR MAIN PROGRAM DTAPE 
(PAGE 1 OF 3) 
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Figure 20 SOURCE AND CROSS REFERENCE LISTINGS FOR MAIN PROGRAM DTAPE 
(PAGE 2 OF 3) 
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Figure 21 TYPICAL DTAPE PROGRAM INPUT DATA SET 



























value to a negative and to control the difference in log exposure between the 
points in the sensitometric calibration data. The difference in log exposure 
between points is computed as 0.15*SENSE and consequently if this difference is 
0.15, the most likely value, the magnitude of SENSE should be set equal to 1.0. 
The subsequent cards in ths data deck contain density values for each point 
which are read into the SWD(I) array. 

The next two major steps in the program involve reading edge identifi- 
cation parameters, the nioraber of points in the edge data and the increment 
between points as well as the edge data itself. These data must also be supplied 
by the user in card deck form. The initial card shown in Figure 21 contains 
two integer identification parameters, the integer NP equal to the number of 
points in the edge data and DX, the increment between points. In the Apollo 
photographic evaluation the identification parameters were set equal to the 
frame number and a sequential number identifying the edge trace made on that 
frame of photography. This card is followed by a sequence of cards containing 
the actual edge data. Provision is made for multiple edge traces in any one 
edge data set. After the last trace a blank card must be inserted to separate 
edge data sets. The user is free to add other edge data sets to the input by 
following this blank card with the new sensitometric calibration data and the 
associated edge data. 

After reading the edge data, the program converts the density values 
to relative exposure by calling function subprogram EXPOS. The converted edge 
data are written onto an output tape and the program recycled to read the next 
edge trace identification parameters and data. This process is repeated until 
a zero is encoimtered for the number of points in an edge trace. This identifies 
the end of the current data set and the program places an end of file on the 
output data tape. It then recycles to read new sensitometric calibration data 
for the next data set. If no data is encountered, the program writes a dummy 
output record as the last file on the output tape and terminates execution. 
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2.0 Function Subprogram EXPOS - This program is part of the Edge Gradient 

Spectral Analysis software package. It supports Program DTAPE in the 
creation of an input data tape for subsequent processing by the main EGSA 
program. The purpose of this function subprogram is to convert a density 
value supplied as input into a relative exposure value using sensitometric 
calibration data. No direct input data must be furnished by the user. The 
input data is supplied through the calling program. No printed output data 
is furnished by this program. 

Figure 22 is a listing of program EXPOS. The function of the major 
blocks are identified by the comments included in the coding. The computing 
path branches depending upon the difference between the input value of density 
to be converted and the density value of the current point in the sensitometric 
data. That point is identified by the integer variable i. If these two 
density values are equal, the exposure for that point is computed and returned 
to the calling program. If the input density value is greater, the value of 
I is increased by one and a comparison between the input density and the 
density value of the next point in the step wedge is made. This process is 
repeated until the density of the point in the sensitometric calibration 
is greater than the input density or until the last point of the data is 
reached. If the latter occurs, the exposure for this last point is computed 
and returned to the calling program. If a point in the sensitometric calibration 
data of greater density than the input value is found the program performs a 
linear interpolation between that point and the previous adjacent point which 
bracket the value of the input density. If the current point in the sensito- 
metric calibration data has a value which exceeds the input density, a check 
is made to determine ^if the density value lies between the current point and 

the adjacent point and if not, the index I is decreased by one and the 

./ 

comparison made again. This process is repeated until the input density value 
is bracketed by two points in the sensitometric calibration. Should the input 
density value be less than the first point in the calibration data (lowest 
density value) the relative exposure is set equal to one and returned to the 
calling program. 
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The relative exposures supplied by this program are 1.0 or greater 
for negative imagery and between 1.0 and 0.0 for positive imagery. The points 
in the sensitometric calibration data must be provided in order of increasing 
density. 

3.0 Main Program EGSA - This program is part of the Edge Gradient 

Spectral Analysis software package. It is one of two main programs in the 
package and is used to compute an Optical Transfer Function (OTF) from edge 
trace data. The complex OTF is expressed as modulus and phase functions. 

The edge data is read from an unformatted data tape created by program DTAPE 
described previously in Section 1.0 of this Appendix. In addition to these 
data, the user must supply card input which identifies various data options 
he elects to receive from the program. Both printed and plotted output data 
are furnished as output options . 

Figure 23 is a listing of program EGSA. The function of the major 
blocks are identified by the comments included in the coding. A macro flow 
diagram for this program was presented as Figure 3 on page 7 of the main text . 

In order to exercise this program, sixteen subprograms are required. 

Four of these programs are user supplied and are discussed in detail in the 
following subsections. These include CONVLV which convolves an input function 
with an 11-point differentiating filter, DXSCAL which determines the increment 
or scale factor (inches/point) for the x-axis of several plots, FORINV which 
performs an inverse Fourier transform and FOURTR which calculates a Fourier 
transform of a real input function providing the complex output function as 
modulus and phase. 

Of the twelve remaining subprograms, nine are used to generate graphic 
output. A brief description of the purpose of each of these programs is 
presented in Section 2.2 of the main text. Since these programs are peculiar 
to the data processing facility, no further discussion is provided in this 
Appendix. Two of the subprograms are IBM supplied FORTRAN IV Library subprograms 
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•• 

•* MAIN PROGRAN - EGSA (EDGE GRADIENT SPECTRAL ANALYSIS) 

• • 

•* THIS PROGRAM COMPUTES AN OPTICAL TRANSFER FUNCTION lOTFI FROM 

*• EDGE TRACE DATA. THE COMPLEX OTF IS EXPRESSED AS MODULUS. ‘TF>, 

•* AND PHASE. *PH>. THE EDGE DATA IS READ FROM AN UNFORMATED TAPE. 

•• DIFFERENTIATED AND FOURIER TRANSFORMED TO COMPUTE THE OTF. INDIVIDUAL 
•• MODULI CAN BE AVERAGED IF MULTIPLE EDGES ARE CONTAINED IN A SINGLE 
•• DATA SET. THE INDIVIDUAL EDGES. THEIR DERIVATIVES (OR LINE SPREADS) 

*• AND RESULTING MODULI ARE PLOTTED IF REQUESTED BY THE USER. OTHER 
•• OUTPUT OPTIONS INCLUDE PLOTS OF THE AVERAGE MODULUS AND ITS 95S 
•• CONFIDENCE BAND AND THE AVERAGE LINE SPREAD AND CONFIDENCE BAND AS 
*• NELL AS PUNCHED CAROS OF THE AVERAGE AND INDIVIDUAL MODULI. 

•• 

•• 

•• *• VARIABLE LIST •* 


•• 
•• 
•« 
•• 
•« 
•« 
•• 
•• 
•• 
•« 
•• 
«* 
«• 
«• 
• • 
•• 
•• 
•• 
•« 
• * 
• • 
•• 
•• 
•« 
•« 
«• 
•« 
•• 
•• 
•• 
•• 
• • 
• • 
• • 
•• 
•« 
•• 
«• 


INPUT VARIABLES 

OX ^INCREMENT BETWEEN EDGE DATA POINTS IN MM. MUST REMAIN 
CONSTANT IN ANY ONE DATA SET. 

E •EDGE TRACE DATA (EXPOSURE VALUES). 

FCO >CUTOFF SPATIAL FREQUENCY USED IN SMOOTHING THE LINE 
SPREAD. 

FINV aSPATlAL FREQUENCY INTERVAL REQUESTED BY THE USER IN 
CYCLES/MM. 

FMAX •MAXIMUM SPATIAL FREQUENCY OF OTF DESIRED BY THE USER. 
EXPRESSED IN CYCLES/MM. 

NED ■OPTIONAL EDGE TRACE IDENTIFICATION VARIABLE. 

VALUE REPRODUCED ON PLOTTED AND PRINTED OUTPUT. 

NIO ■OPTIONAL EDGE TRACE IDENTIFICATION VARIABLE. 

NP ■NUMBER OF POINTS IN EDGE TRACE DATA) < OR ■ 900. 

OPTION ■ALPHA NUMERIC CHARACTERS IDENTIFYING USER SELECTED 
OUTPUT OPTIONS. 

MAJOR INTERNAL VARIABLES 

ALS bAVERAGE line SPREAD FUNCTION. 

ALSM ■AVERAGE LINE SPREAD CONFIDENCE BAND I LOWER BOUND) . 

ALSP ■AVERAGE LINE SPREAD CONFIDENCE BAND (UPPER BOUND). 

ATF ■AVERAGE MODULUS OF THE OTF (OR MTF). 

ATFM >AVERAGE MTF CONFIDENCE BAND (LOWER BOUND). 

ATFP bAVERAGE MTF CONFIDENCE BAND (UPPER BOUND). 

DNU bACTUAL SPAtlAL FREQUENCY INTERVAL OF OTF IN CYCLES/MH. 

KNT bCOUNTER equal TO THE NUMBER OF EDGE TRACES IN A DATA SET. 
NCF bNUMBER of points IN DIFFERENTIATED EDGE TRACE. 

NPL bNUNBER of INDIVIDUAL FUNCTION PLOTS REQUESTED BY USER. 

NPLS bNUMBER OF POINTS IN THE STORED OR UNSMOOTHED LINE SPREAD. 
NPTF oNUNBER OF POINTS IN THE OTF. (bI^NUMBER OF HARMONICS) 

NSL bNUMBER OF POINTS IN THE SMOOTHED LINE SPREAD FUNCTION. 

PH bPHASE of THE OTF. ON. 

SL bLINE spread function (SMOOTHED VALUE IF OPTION REQUESTED). 
SSL ■■STORED* OR UNSMOOTHED LINE SPREAD FUNCTION. 

TF bMOOULUS OP THE OTF. 
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ISN 0002 


ISN 0003 
ISN 0004 
ISN 0005 
ISN 0006 
ISN OOOT 
ISN 0008 
ISN 0009 


ISN 0010 
ISN 0011 
ISN 0012 


ISN 0013 
ISN 0014 


ISN 0015 
ISN 0016 


ISN OOIT 
ISN 0018 
ISN 0019 


c 

** 


c 


AUTHORS: R. E. KINZLV 

c 


M. J. MAZUROMSKI 

c 


OATES JULY. 1966 

c 

•• 

REVISIONS: AUGUST. 1973 

c 




c 

OINENSION XP<900I.ATF(900ltETFI900t,SLI900)>EI900I.CFI9CO),TF(900) 
*iATFMI900ltATFP(900ltl>HI900I.DUMIU,XPI(900>.XLII90C>.OPTION(T ), 
•ALSI900I.ELSI900lfALSPI900I.ALSN(900ltXL(900)«SSLI900l. 

•XS( 9001 .OPT( 71 ,XLE(900I 
LOGICAL LTV(TI/7*. FALSE./ 

DATA QPT/4H E.4H Lt4H T.4H SLi4N AL>4N AT.4H PT/ 

2000 FQRNAT(2F10.0tTA4.FI0.0l 

3000 FaftNAT(//10X,2110 / UOXtlOFIO.3) I 

3010 FORNATIIIOI 

3020 F0RHATIF5.3>llF6.3i 

3030 F0RNAT<lH1.7Xf • > THE Fft£Q(£NCV INTERVAL OF THE OTF HAS BEEN CHA 

*NGEO FR0N'.F6.2.* CVCLES/NN TO*t F6.2, • CVCLES/NH, • / 13X. 'THE 
*EOGE OATA SET IDENTIFICATION PARAMETERS ARES NED* '.HO.'. NID* •. 

*110 I 

3032 FORNATiaXt • USER REQUEST FOR AN AVERAGE LINE SPREAD HAS BEEN 

•OELEATEO.' I 

3034 FORMAT! 8Xt • > USER REQUEST FOR AN AVERAGE NTF HAS BEEN OELEATE 

♦O.M 

CALL PLOTERU.l.lOl 
C 

c «**•««•••••*««•*«••••««*«••«*•**•*•••••••••««**«**•*••••••»**••••••*»•*•*••** 

c * • 

C • READ USER INPUT OATA AND SPECIFIED OPTIONS * 

C • ■ * 

C a**************************************************************** •*••** *»•**• 

C 

10 REA0IS.200O.ENO>5S0IFINV.FMAX,IOPTIONIIItI«ltT>,FC0 
NPL*0 


C 

C 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


•*••«••••*•*« p*************** •*••*«••••••••••••••••*•***••***•*•• •*••**•*•«*• 

* * 

* INITIALIZE LOGICAL ARRAY 'LTV* DESCRIBING OUTPUT OPTIONS * 

• * 


DO 16 I*lf7 
16 LTVIIi*. false. 

• * 

• SET 'LTY* array ELEMENTS TO 'TRUE* FOR OPTIONS REQUESTED « 

* * 
••*•«•*•••••••*•*«••••«••••**•*•••**•*•*•*** •*•••••* *«••** 


00 IB 1*1.7 
00 18 J*lf7 

18 IF ( OPTION! 1 1. EQ. OPT UlILTVfJI*. TRUE. 
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ISN 0021 
I$N 0022 
ISN 0023 
ISN 0024 
ISN 0025 
ISN 0026 
ISN 0027 


ISN 0020 
ISN 0029 

• ■ i ■ 


ISN 0031 
ISN 0032 
ISN 0033 


ISN 0034 
ISN 0035 
ISN 0036 
ISN 0037 
ISN 0038 


ISN 0040 


Figure 23 
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****«***«***«*«6*****««**«««« «••••*•••••*•••••«*••••••••••••••••••••••« «**••• 

* ■ ■ • 

• B6GIN PROCESSING CURRENT EDGE OAT* SET • 

* . ■ • 

NPTPbFNAX/F1NV«I . 5 
0F«FH*X/5. 

20 KNTbO ' 

CALL CLE*R(ATFmt*TFI900ll 
CALL CLE*RIETFUIfETF(900ll 
CALL CLE*R(*LSat<*LSI900n 
CALL CLEAR(ELSIlltELS(900l> 


• «••**•••«••**•*••• 6«««* •«••••««•«••••*• *«*««*»»**P****P***«*»*** •*•••• •••«•• 

• • 

* READ EDGE DATA PARAMETERS • 

•*•«**••••••• *«»***6*« 

50 REA0llfERR«550tEN0a400INE0tNI0tNPt0X 
IFINP.EO.OI GO TO 400 


**«••*•••*«••••«•«• *««*««««««**«»**64****»««*6P**6P**6**P*M**********» •••••• 

• ■ • 

« READ EDGE DATA FROM TAPE » 

••**•**«*«••«*«•* «•****•«••* *•*••• 

NP»NP*12 

REAOIl t ERR«SSO tENOaSSOI I E (I I tI-13 tNP I 
KNT«KNT*1 

••••••••••••A** •«*«••«*••• •«•*••••••••**••*• *•***• 

• • 

• ADJUST EDGE DATA FOR NUMBER OF POINTS TO BE LOST IN DIFFERENTIATION • 

• • 


DO 55 l-ltl2 
EIII«E<13I 
55 E(NP*II«EINP| 

NP*NP^12 

IFILTVtlll CALL SCALEIE.NP.2.tOUNNy.OUMMY,II 


*6«** •*•*••••**•*•*• *****«**»«*i 

• DIFFERENTIATE EDGE DATA 

•*••••••••• »«««««««6««*««*««**««««*«*M*6*««**«*****«*6»******«** *•••••*••** 

CALL C0NVLVIEtNPtDUHtIfCFfNCF»3l 
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ISN 0050 
ISN 0052 
ISN 0053 
ISN 0054 
ISN 0055 
ISN 0056 


ISN 0057 
ISN 0059 
ISN 0060 


ISN 0061 
ISN 0062 
ISN 0063 
ISN 0065 
ISN 0066 
ISN 0067 


Figure 23 


c **•***•••••••••••••**••*••**••••*•*•••*•••••*•*•«***••*•**»•**•***•**•***•»» 

c • 

C * C0MPU7E THE NUMBER OF POINTS NEEDED IN LINE SPREAD TO OBTAIN DESIRED 
C * FREQUENCY INTERVALt 'FINV't FOR OTF 
C • 

c •*••*«***••***•••*•••••••**«***•••**••*••••«•••«*••*••**•*•«•*•*»••••«****•* 

c 


c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


NPLScl./(FINV*DXI«0.5 

* 

• SET UP X-AXIS INTERVALS AND ARRAYS FOR INDIVIDUAL EDGES AND LINE 

• SPREADS AND AVERAGE LINE SPREAD IF REQUESTED BY USER 

• 

• A********************************** ••**•*•** •*«•••*•••*••*•*»•*•***••*«*••« 

IF(.NaT.(LTV<SI.OR.LTV(2>.OR.LTV<UI> CO TO 70 
IFtXNT.EQ.II CALL OXSCALI2.5tNPLSiDXtDXTItOXl I 
CALL OXSCAL(2.5,NP«DX.OXTE.OXEI 
DO 60 I'lfNPLS 
XL«II=( 1-I-NPLS/2190X 

60 XLE(II>(I-1I*0XE 
IFIKNT.NE.il 60 TO 62 
DO 61 I‘I.NPLS 

61 XLHII«(I-1I*0XI 

62 CALL SCALE(XLtNPLS.6..XN.DItli 
70 CALL CLEAR(SL(Ii.SL(NPLSII 

CALL CLEARISSLIlltSSLINPLSn 


• ' * 
« TEST FOR OECREASING EDGE DATA AND ADJUST IF NECESSARY « 

• • 
**«*•«•«*««*«•*«•«« «**•**••«* •*»•«•*«*••**•••**•«••••*•*••*•**•»••••••* *••*•• 


IFIEINPI-Elll.GT.O.I GO TO 72 
DO 71 I«1.NCF 
71 CFlIle-CFlII 

• LOCATE THE NAXINUN VALUE OF THE DERIVATIVE OF THE EDGE DATA * 

• « 

•*«••*••*•••*•*••«• a******************** ••••*••«•••••••*••*•*•«***•** ****«••• 


72 BGLaO.O 

DO 80 I'ltNCF 
IFICFIII.LT.BGLI GO TO 80 
BGL>CF(II 
LBIG«I 
80 CONTINUE 

a*************************************** *•**•*••••••****••••****•••**** aaaaa 

a 

a LOCATE AND SMOOTH THE NEGATIVE VALUES IN THE TAILS OF THE 
a DERIVATIVE. IF ANY 

a 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa aaaaa 
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ISN 0068 
ISN 0069 
ISN 0070 
ISN 0072 
ISN 0076 
ISN 0075 
ISN 0077 
ISN 0078 
ISN 0079 
ISN 0080 
ISN 0082 
ISN 0086 
ISN 0085 
ISN 0087 


ISN 0088 
ISN 0089 
ISN 0090 
ISN 0091 
iSN 0092 
ISN 0093 
ISN 0096 
ISN 0096 
ISN 0097 
ISN 0098 
ISN 0099 
ISN OiOt 
ISN 0102 
ISN 0103 
iSN 0106 
iSN 0106 
ISN 0107 
ISN 0109 
ISN 0110 
ISN 0111 
ISN 0112 
ISN 0113 
ISN 0116 
ISN 0115 
ISN 0117 
ISN oils 
ISN 0120 
ISN 0l21 
iSN 0122 


I7EST<:LBIG 
90 ITEST=ITEST*1 

IFIITEST.GT.NCF-II GO TO 101 
IFICFCITESTI.GE.O.i GO TO 90 

CF(ITEST»««CFaTEST-U*CFIlTESTJ*CF(ITEST*ll»/3. 

IF(CF(1TEST).GE.0.> GO TO 90 
101 IHI=ITEST-1 
ITEST^LBIG 
105 ITEST-ITEST-1 

IFIITEST.lt. 21 GO TO 110 
IF(CF(ITESTI.GE.O.) GO TO 105 

CFI1TEST»»«CF«1TEST-1I*CF(ITESTJ*CFUTEST»1)I/3. 

IFICFIITESTI.GE.O.) GO TO I OS 
110 1L0«ITEST*1 
C 

C ••**•«•••• •••••**••••«*«*••••«• «•••••*••«***•» •• •**••••••• 

C • ' ' * 

C • ADJUST THE NUMBER OF POINTS IN THE DERIVATIVE TO INCLUDE THE PORTION • 

C • BETWEEN 'ILO* AND* IHI* AND EQUATE THE LINE SPREADS TO THIS PORTION • 

C • • 

c 

c 

NCF>1HI-IL061 

112 KP«INPLS*U/2 
KLO^LBIG-ILO*-! 

KHI>1HI-LBIG 

1B«KP-KL0 

1E«IL0-1 

IFIIKLO.GT.KPl.OR.IKHI.GT.NPLS-KPII GO TO 116 
00 113 I«lfNCF 

113 SLI1B«I)*CFIIE*11 
GO TO 116 

116 IFIKLO.GT.KHll GO TO 115 
NPLS-29KH1 

FlBl.O/(FLaAT(NPLSi«OXI 
WtlTE«6t3030> P1NV^FI,NED>N10 
IFILTVI5II WRITE «6|3032I 
LTVI5I>. FALSE. 

IFiLTVI6ll WRITE I6i3036l 
LTV(6I«.FALSE. 

FINV«F1 
GO TO 112 
115 NPLS«2«KL0 

FI«1.0/IFL0ATINPLSi*0X) 

WRITEI6. 30301 FINV.fi, NED.NIO 
IF(LTV(5II WRITE 16,30321 
LTV(5t>. false. 

IF(LTVI6II write 16,30361 
LTVI6)>. FALSE. 

FINV«FI 
GO TO 112 
C 

C o*************************************************************************** 
C • • 

C * NORMALIZE LINE SPREAD AND EQUATE TO "STOREO* LINE SPREAD * 

C * • 

C •••••***«***•*«****««»**•*•*•****•••*«•••*••*•*•*•*•*•*••******•••*•••****«*« 
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ISN 0123 
ISN 0124 
ISN 0125 
ISN 0126 
ISN 0127 
ISN 0128 


ISN 0129 
ISN 0131 
ISN 0132 
1$N 0133 
ISN 0134 
ISN 0136 


ISN 0138 
iSN 0139 


ISN 0141 . 


ISN 0142 
ISN 0143 
ISN 0144 
ISN 0145 
ISN 0146 
ISN 0148 
ISN 0149 
ISN 0150 
ISN 0151 
ISN 0153 
ISN 0154 
ISN 0155 


Figure 23 


c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


116 SUM'0.0 

00 117 I'ltNPLS 

117 SUNaSUM«SL(ll 
DO 118 I«1>NPLS 
SL(II«SL(II/SUN 

118 SSL(II>SL(II 

• •«•**•••**«•••«*•*•«••••••*•**••»*•••••**•••*«••*•••***••*•••*•**•••*•***••* 

• « 

* COMPUTE STATISTICS FOR THE AVERAGE LINE SPREAD « 

• ' » 


IF(.NOT.LTVISI) GO TO 130 
DO 120 I>1»NPLS 
ALSIIi>ALSIll«SLIll 
120 ELSm°ELS(I)*SLIII*SL(I> 

130 IFILTV(2I)CALL SCALE(SSLiNPLSt2. tDUHNrtOUHMV.l I 
lF(.mT.ILTVI3).0R.LTVI4).0R.LTV(6tll GO TO 210 

A****************** •«••«•• «»••«*«•*••««•*«••«*•*«»•*••••••***••**•••*»• *«*«** 

• • 

• FOURIER TRANSFORM LINE SPREAD TO OBTAIN THE COMPLEX OPTICAL TRANSFER * 

• FUNCTION * 

• * 
•A********************** •••••**•**•••»••**•*••*«•• •••*••••••*•••**•*••* ***••* 

CALL FOURTRISLtNPLSfOXtTF.PHfDIRI.NHI 
IFI.N0T.LTVI4II GO TO 160 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 

A A 

A COMPUTE THE NUMBER OF HARMONICS FOR THE SMOOTHED LINE SPREAD A 

A A 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 


NHNaFCO/FINV 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 

A A 

A APPLV INVERSE TRANSFORM TO OBTAIN THE SMOOTHED LINE SPREAD a 

A * 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 

CALL FORINVITFtPHtNHM.DNUtSLtNSL.DXSLStll 
DO 140 Ia1,NSL 
140 SLIIIaABSISLIIII 

CALL SCALEISL.NSL. 2. .DUMMY, DUMMY, 11 
IFIKNT.NE.il GO TO 180 
CALL OXSCAL(2.5,NSL,DXSLS,DXSI,OXSI 
DO 150 I°1,NSL 
150 XS(I|a(I- 1IA0XS 
160 IFIKNT.NE.ll GO TO 180 
DO 170 1=1,NPTF 
XPCIl-ll-IIAONU/OF 
170 XPUI)a2.5*XP(1)a0.5 
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ISN 0156 
ISN 0157 
ISN 0158 
ISN 0159 
ISN 0160 


ISN 0161 


ISN 0163 
ISN 0166 


ISN 0165 
ISN 0167 
ISN 0168 


ISN 0169 
ISN 0170 
ISN 0171 
ISN 0173 


ISN 0175 
ISN 0176 


Figure 23 


V*************************************** ••***••••**••*•*••«*•••••*• ••**••**• 
* 

* COMPUTE SUMS RELATED TO THE AVERAGE MTF AND CONFIDENCE BAND 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


180 TFl^TFIll 

DO 200 I«lfNFTF 
TFI 1I«TF(I1/TF1 
ATF(l|-ATFm*TF(II 
200 ETF(II«ETF(II«TFIII*«2 

v*********** 

« • 

• PUNCH INOIVIOUAL MTF IF REQUESTED BV USER • 

• . • 


1F(LTV(7II HRlTE(6,3000)NEO.NIOtfrF(IltI*l>NPrFI 
• SCALE INDIVIDUAL NTF FOR PLOTTING • 

6 - V 


DO 205 I«1»NPTF 
205 TFItlaTFIlMI.O 

• PLOT INOIVIOUAL FUNCTIONS AS SPECIFIED BV USER 


210 IFIMQOIKNT.lSI^E.ll GO TO 216 
FNEOaNEO 
NOSaNEO 

aaaaaaaaaaa«a«aaaaaaaaaaaaaaaaaa*aaaaaaaaaaaaaaaaaaaaaaaaaaaaaa*aaaa**aaaaaaa 
a • 

a ESTABLISH THE NUMBER OF INDIVIDUAL FUNCTION PLOTS REOESTED a 

a a 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 


NPLaO 

00 212 Ial,6 
212 IFILTVUII NPLaWPL*! 

IFINPL.EQ.OI GO TO 50 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 
a ' ' a 

a ORAM AXIS FOR EACH TYPE OF INDIVIDUAL FUNCTION PLOT REQUESTED a 

a a 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 


CALL PLOTIO.fO.tOI 
XPLa-B.5 
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ISN 0177 
ISN 0178 
ISN 0179 
ISN 0180 
ISN 0181 
ISN 0182 
ISN 0183 
ISN 0189 
ISN 0185 
ISN 0186 
ISN 0187 
ISN 0188 
ISN 0190 
ISN 0191 
ISN 0192 
ISN 0199 
ISN 0195 
ISN 0196 
ISN 0197 
ISN 0198 
ISN 0199 
ISN 0201 
ISN 0203 
ISN 0205 
ISN 0207 
ISN 0208 
ISN 0209 
ISN 0210 
ISN 0211 
ISN 0212 
ISN 0213 
ISN 0219 
ISN 0216 
ISN 0217 
ISN 0218 
ISN 0219 
ISN 0220 
ISN 0222 
ISN 0223 
ISN 0229 
ISN 0225 
ISN 0226 
ISN 0228 
ISN 0230 
ISN 0232 
ISN 0233 
ISN 0239 
ISN 0235 
ISN 0237 


ISN 0239 
ISN 0291 


DO 219 I^l.NPL 
XPL=XPL*8.5 

CALL MGRI0(O.S.XPL*5.6,1O..2.5.2..2.5i1> 

CALL NGRID(0.S.XPL*3.0tl0.,2.5.2..2.5tII 
CALL MGRID(0.S.XPL»0.9.1O.t2.5t2..2.5,l> 

219 CALL NUNBERI0.1SiXPL«9.S.0.19.PNED»-90..-ll 
XR°0.0 

CALL PL0TI8.5i5.6.-3i 
XOR 1=0.0 
YOR-0.0 
GO TO 218 

216 IFINPL.EQ.OI GO TO 50 
YOR«0.0 
XDRo-2.0 

IF(N00(KNTt5I.NE.ll GO TO 218 

X0R«8.0 

TOR >-2.6 

218 ULL PL0T(X0RtY0R*XR|-3l 
XR«0.0 
NPLT>NPL 

IF I LTV! in GO TO 230 

221 1FILTVI2I) 60 TO 231 

222 IF<LTVI3n GO TO 232 

223 IF(LTV(9I) 60 TO 233 

230 (••LT-NPLT-1 

CALL LINEIEtXLEtNPdl 

CALL NUNBER(0.8t2.9.0.0T«OXTE>-90.t3l 

KC«1 

GO TO 290 

231 l#*LT«NPLT-l 

CALL LINEISSLtXLItNPLStll 

IFIKNT.EQ.ll CALL NUMBERt0.a,2.6.0.0T>DXTl ,-90.«3) 

KC«2 

GO TO 290 

232 NPLT«NPLT-1 

CALL LINEITFtXPltNPTFtll 

IFlKNT.EQ.il CALL NUNBER(1.8t0.T5t.0T.FNAX,-90.tll 
KC>3 

GO TO 290 

233 NPLT*NPLT-1 

CALL LINEISL.XS.NSL.ll 

lFIKNT.EQ.il CALL NUN8ERI0.8.2.9t .0T.DXSI.-9O..3I 
IFIKNT.EQ.ll CALL NUHBER(1.0»2.9.0.0T.FCO.-90.tll 
290 IFINPLT.EQ.OI CO TO 90 
CALL PL0TI0..8.5.-3I 
XRAXR-B.S 

GO TO l221.222f223l.KC 
900 IFIKNT.EQ.OI GO TO 10 

1F(.N0T.ILTVI5I.0R.LTVI6III GO TO 20 
C 

c *•«•••••••««••««•«•••«••••••••••••••••••••••«•«•••*•••••*••*•*•*****••*•*•*•* 

c • • 

C • CALCULATE AVERAGE NTF ANO 95S CONFIDENCE BAM) • 

C • . • 

c •••«*****«*«**«««*««««*««»«*«*****«***««*«*** 6 ** 9 ******»****«******«*** 

C 

IF) .N0T.LTVI6II GO TO 960 
00 950 lol.NPTF 
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tSN 0242 
ISN 0243 
tSN 0244 
ISN 0245 
ISN 0246 
ISN 0247 
ISN 0248 


ISN 0250 
ISN 0251 
ISN 0252 
iSN 0253 
ISN 0254 
ISN 0255 
ISN 0256 
ISN 0257 
ISN 0258 
ISN 0259 
ISN 0261 
ISN 0262 
ISN 0263 
ISN 02M 
ISN 0265 
ISN 0266 
ISN 0267 
ISN 0268 
ISN 0269 


ISN 0271 
ISN 0273 
ISN 0275 


ISN 0277 
ISN 0278 
ISN 0279 
ISN 0280 
ISN 0281 
ISN 0282 
ISN 0283 
ISN 0284 


Figure 23 


ATFIII«ATF(Ii/KNT 

DEV«ANAXlllETFIII-«NT«ATF(II*«2ltO.OI 
0EV>>2*SQRT(DEV/(KNT-1I I 
OEV>OEV/SORT(FLOAT(KNT) I 
ATFNIIi>ATFIII-OEV 
450 A7FFm«ATFIII«0EV 
460 IFI.NOT.LTVI5I) GO TO 490 

C ••4464446646446M496M44M6666*********************************************** 

c • • 

C • CALCULATE AVERAGE LINE SPREAD AND 95t CONFIDENCE INTERVAL • 

C • * 

c **«*«*4666*6«66«««*****«*6*6«6*«»**6«»*6*66*«************************«* *••••* 


00 480 I>ltNPL$ 

ALSIII>ALSIII/KNT 

0EV>ANAX1IELS 1 1 1 HUITVALSII l««2«0.0l 
0EV«2*SQRT(0EV/( KNT-1 I I 
OEV«OEV/SORT(FLOATIKNTI I 
ALSP(I>«ALSIIl40EV 
480 AL$N(II-ALSIII-OEV 
BGL-O.O 

00 482 I*1«NPLS 
IFIALSPI1I.lt. 86LI GO TO 482 
B6L>ALSP|ll 
LBIG«l 

482 CONTINUE 

CALL SCALEIALSP.NPLStS.tSLNtOL.ll 
SP«ALSP(LB1GI/B6L 
DO 484 I«1,NPLS 
ALSIII-AUUIPSF 
484 ALSNIII«AL$HUI*SF 

1FI.N0T.LTV(6II GO TO 300 
C 

c 666666666666666666666966M6A66666666666669999P******************************* 

c • • 

c • PUNCH AND PRINT AVERAGE NTP IF REQUESTED * 

C • • 

C 66666666666669M66666666666666666666P66P6P666666P6PP6PPPP6************* 

c 

490 IFILTVITII HRITE (6.3000IN0S.KNT.I ATFI li.I>l.NPTF I 
IFILTVmi HRITEI7.3010INDS 
IFILTVITII HRITE 17,30201 I ATFIll.I-l.NPTF) 

C 

C •66»»»66666666»66»6»66—69» »66— 6 6 »9>6»6 6666* »—>»»«»»*»»»**»»»»»»»*»«* **«»** 

C 9 • 

C • PLOT AVERAGE NTF AND CONFIOBICE INTERVAL IF REQUESTED • 

C • • 

C •66666666966444466666666M66996666666666******************************* 


DO 495 I«I,NPTF 
ATFII)-5.0*ATFIII 
ATFPIII°S.O*ATFPIII 
495 ATFNII|s5.09ATFNII| 

CALL PLOTIO.,0.,01 
CALL PL0TI2.5,6.5,-3I 

CALL AXISI0.,0.,2HNU,-2,5.,-90.,O.,0FI 
CALL AX IS 1 0 • f 0* f 3HTAU f 3 f 5* f 0*t0*t0*Z) 
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ISN 02BS 
ISN 0286 
ISN 0287 
ISN 0288 
ISN 0289 
ISN 0290 
ISN 0291 
ISN 0292 


ISN 0294 
ISN 0295 
ISN 0296 
ISN 0297 
ISN 0298 
ISN 0299 
ISN 0300 
ISN 0301 
ISN 0302 
ISN 0303 
ISN 0304 
ISN 0305 
ISN 0306 


ISN 0307 
ISN 0308 
ISN 0309 


C 

C 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


CALL HGRI0(0.i-5.t5.i5.tI.tI>tU 
CALL LINE(A7FP.XP.NP7F«1> 

CALL LINEIATFiXP.NPTF.ll 
CALL LINE! ATFN.XPtNPTFtll 

CALL SYHBOL(-I.S.-1.2.0.14.22HAVERAGE OF SAMPLES,-90.,22 I 
CALL NUHBERI-1.5.-2.S2.0.14,FL0ATIKNT>.-90.t-ll 
CALL NUMBER (5. 5. -2. 16,0. 14. FNED. -90.. -1) 

IF(.N0T.L7V(5I I GO 70 20 

• • 

• PLOT AVERAGE LINE SPREAD AND CONFIDENCE INTERVAL IF REQUESTED * 

• • 


500 CALL PL0T(0.,0..0l 
CALL PLOT(2.5,7. ,-3l 
CALL AXIS<0.,0.,lHX.-1.6..-90..XM,0Ii 
CALL AXIS(0.,0..4HLIXI.4,5., O..SLM,OL> 

CALL PL0T(0.,-6.,-3l 
CALL LINElALS.XL.NPLS.ll 
CALL LINEIALSP.XL.NPLS.li 
CALL LINEtALSMtXL.NPLS.il 
CALL PL0T(0..6..-3I 

ULL SVMB0L(-l.S.-l.T.0.t4.22HAVERA6E OF SAMPLES. -90.. 22 > 

CALL NUMBER (-1.5 .-3.02. 0.14,FL0ATIKNT I .-90. .-1 I 
CALL NUMBERI5.S.-2.66tO.14.FNE0.-90..-U 
GO TO 20 


• ' * 

• TERMINATE EXECUTION * 

• * 
A****************** •••«•«• •*•*«•••••*•••*•••«••••••«•*•*•*«*«****** ********** 

550 CALL EFPLOT 
STOP 
B40 
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L 1 S T I N C«*** 


:;b<: r rOjR iTR a N ,/;lC ROSS REFERENCE 

SYMBOL INTERNAL, statement NUMBERS; . '-1) :V: ■ X T: ' * 


E 

0002 . 

0032 

0C35 

0035j;.^0036 0036 -.'0038 ui0040 

0057 

COST 

0208 








1 

0013 

0013 

oou 

0015.';;u016 ..0017. 0019 

0032 

0032 

0032 

0034 

0035 

0036 

0047 

00^6 

0048 0049 

0049 

0052 


0053 

, 0053 

0059 

^ 0060.7:0060 L :0062 0063 

0065 

0066 

0096 

0097 

0097 

0124 

0125 

0126 

0127 0127 

Ul26 

0128 


0131;; 

: 0132 

0132 

•0132 ^ 0133 , 0133 0133 

-0133 

0143 

0144 

0144 

0149 

0150 

0150 

0153 

0154 0154 

0155 

C155 


0157 

0158 

0158. 

0159 j ,0159; :0159..j 0160 c 

:0160.. 

0160 

0161 

^^^0161- 

^0161r. 

..0163>oc0l64.'t 

iC164o 

c0l70 ^0171 

0177 

C241 


02A2 

0242 

0243. 

0243 0246 0246 0247 

0247 

0250 

0251 

0251 

0252 

0252 

0255 

0255 

0256 o 0256 

0258 

0259 


0261 

0262 

0266 

0267,1 0267 0268 0268 ' 

0271 

0271 

,0271 

0275 -0275 • 

'0275.x 

•0277 

0278 

0278* 0279 

0279 

0280 


0280 












c 



J 

0016 

0019 

0019 

- " ■ . i'.- .. • : j-,-; 

' 

' ^'•'5 




00 -r- 

. *•::• , 



CF 

0002 

0040 

0060 

0060 0063 0065 0072 

0074 

0074 

0074 

0074 

0075 

0062 

0084 

0084 

0084 0084 

0085 

0097 

DF 

0022 

0154 

0283 













01 

0054. 

0296 














OL 

0264 

0297 


. r ...... 












OX 

0028 

0041 

0044 

0046 0048 . 0102 0113 '. 

.0138 











FI 

0102 

: 0103 

0110 

0113 .0114. 0121; ; 












IB 

0092 

0097. 














IE 
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ORIRAN CROSS REF6RENCE LISTIN &•**•* 
SYHBOL INTERNAL STATEMENT NUMBERS 
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and should be available to most users. The last subprogram, CLEAR, fills a 
specified block of storage with zeros. If necessary, an equivalent operation 
can be accomplished by a sequence of do loops in Program EGSA. 

The program initially reads user supplied input data on punched cards 
The following variables are specified for each data set to be processed: 

FINV - spatial frequency interval desired for output OTF data in 
cycles/mm 

FMAX - maximum spatial frequency of output OTF data in cycles/mm 
OPTION (7) - alpha-numeric array identifying the types of output 
data desired by user 

FCO - cutoff spatial frequency in cycles/mm used in smoothing the 
line spread 

Seven output options are available to the user. These are specified by 
employing one or more of a set of alpha-numeric characters on the punched 
card input. These characters must be right hand justified in the appropriate 
field and may be listed in any order. The following characters are allowed: 

E - specifies plots of the individual input edge traces of the data 
set 

L - specifies plots of the individual line spread data resulting 
from differentiation of the edge data 
T - specifies plots of the individual MTFs* resulting from Fourier 
transformation of the line spread functions 
SL - specifies plots of smooth line spread functions 
AL - specifies a plot of the average line spread and its associated 
confidence bands 

AT - requests a plot of the average MTF and its confidence 
bands 

PT - requests the values of the individual and average MTFs 

* The Modulation Transfer Function (MTF) is simply the modulus of the complex 
OTF. 
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on punched cards and printed values of the average MTF. The 
values of the average MTF are not punched or printed unless the 
option for a plot is also specified. 

The program sets elements in. a logical array, LTV, to true for the options 
requested if they are properly specified on the input data card. Figure 24 
contains an example of this data card where the user has selected plots of the 
individual MTFs and the average MTF. This input card also specifies the 
spatial frequency interval as 5 cycles/mm and a maximum spatial frequency of 
150 cycles/mm. The cutoff spatial frequency variable is unspecified or equal 
to zero, indicating that the user has not elected to smooth the line spreads 
consistent with the deletion of the SL option on the card. The value for 
the frequency cutoff is, in fact, ignored if this option is not specified. 

The next major step in the program is the initialization of parameters 
and clearing of arrays prior to beginning the actual processing of an edge 
data set. Identification and data: records are then read from an input data 
tape created using main program DTAPE described previously in Section 1.0. 

The edge data are then adjusted for the number of points that will be lost 
due to differentiation. The differentiation is accomplished by subroutine 
CONVLV described in the next section. Although the current version contains 
an 11-point, internally specified differentiating filter the dimensions allow 
a maximum of a 25-point filter. The adjustment in the edge data is commensur- 
ate with the maximum size filter. 

After differentiation the next major step in the execution involves 
the set-up of x-axis arrays for the output plots of several of the individual 
functions. This is followed by manipulation of the differentiated edge data 
including a non-stationary 3-point smoothing of negative values in its tails. 
Only the non-negative portion of the edge data is retained and equated to the 
smoothed and unsmoothed line spread functions. The line spread function is then 
Fourier transformed by calling FOURTR to obtain the modulus and phase functions 
of the OTF. If the user has requested a smooth line spread, this is computed by 
truncation of the OTF and subsequent inverse Fourier transformation. The 
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Figure 24 TYPICAL EGSA PROGRAM INPUT DATA CARD 

















inverse transformation is accomplished by subroutine FORINV. This calculation 
is primarily a cosmetic operation and produces no new information about the 
optical system performance that is not contained in the unsmoothed line spread 
or the complex OTF. 

The next steps in the program involve the production of output data 
as requested by the user. These include punched cards of the individual MTF 
and plots of the individual edges, line spreads, smoothed line spreads or MTF. 

The program' then recycles to read the next edge trace in the current edge data 
set and repeats the processing cycle. If no additional edge data exists, the 
calculation and plotting of the average line spread and a 95% confidence band 
and the average MTF and its confidence band are accomplished if these options 
were requested by the user. The program recycles to process the next edge data 
set, if any, in the same manner. If no additional data sets are identified 
the program terminates. 

Examples of the plots of individual function data generated by the 
program are presented in Figure 25. Part (a) shows individual line spread 
functions which resulted during an analysis of the reproduction of Apollo 15 
photography and Part (b) shows individual MTF data resulting during a similar 
analysis of Apollo 17 photography. The original plots fill an 8 x 11 1/2 
inch page and were reduced for reproduction here. Similar plots are generated 
for the individual edges and smoothed line spreads if requested. A number is 
placed at the bottom of each page which corresponds to the identification number 
NED specified by the user.. In the examples shown this is the number of the 
frame from which the edge data was acquired. The top left hand plot in both 
cases also contains a number specifying the total length of the x-axis. In 
the case of edges, line spreads or smoothed line spreads, this length is 
expressed in millimeters and for MTFs as spatial frequency in cycles/millimeters. 
Figure 26 contains an example of an average MTF plot generated from the 
individual MTFs in Figure 25(b) as part of the evaluation of Apollo 17 photo- 
graphy. The upper and lower curves define the 95% confidence band. The number 
placed at the top of a plot corresponds to the user specified variable NED and 
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Figure 26 EXAMPLE OF AN AVERAGE MTF PLOT FROM THE EGSA PROGRAM 
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in this case is equal to the corresponding frame number. The number of 
individual edge traces which were used to compute the average MTF is indicated 
on the bottom of the plot. A plot with a similar format is generated for the 
average line spread function if that output option is requested. 


4.0 Subroutine CONVLV - This program is part of the Edge Gradient Spectral 

Analysis software package. It supports Program EGSA in the calculation 
of an OTF from edge trace data. The purpose of this subroutine is to differen- 
tiate the edge trace data using an internally specified differentiation filter. 
No direct input data must be furnished by the user. The input data is supplied 
through the calling program. No printed output data is furnished by this 
program. 

Figure 27 is a listing of subroutine CONVLV. The fiinction of the 
major blocks are identified by the comments included in the coding. This 
subroutine has capabilities that are not exploited by the EGSA main program. 

It can perform the convolution between two input functions or an input 
function and an internally stored filter. The internal filter options include 
ah 11-point smoothing filter and ein 11-point differentiating filter. A maximum 
of 25-point filters are permitted by the dimension statement and can be inserted 
by changing the appropriate data statement. The options available in the 
convolving operation are controlled by the variable KCNTRL as follows: 

1 - specifies the convolution between two input functions 

2 - causes the primary input function to be smoothed by the internal 

filter 

3 - specifies differentiation of the primary input function by the 

appropriate internal filter 

4 - indicates smoothing of the primary input function followed by 

differentiation. 

The computing task branches depending upon the value of the variable KCNTRL. 

When the specified option involves the convolution between two input functions 
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SMOOTHING FILTER DEFINITION lOATA STATMENTl 
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Figure 27 SOURCE AND CROSS REFERENCE LISTINGS FOR SUBROUTINE CONVLV 
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SMOOTH INPUT FUNCTION “FUNl- USING THE INTERNAL FILTER -SHF" 



z o 
^ o 


(/> </) 
uo O 
O -j 


• z 

II 3 
> IL. 

i -1 o 


♦ o 

O 

i/> — 

o z — 

-J D 4 - 
vO II. • 

O z U CJ 

-i ♦ D Uj 
» uo u. II • 

1 ^ i/> <_> ^ 

or II o II »» or 

O h- 

Z I i-« z z 

H-* (7“ Lj -» *» Ui qC 

K 4 - Z Z Z ^ D 

Z U 3 S «- I- O 

OQi^U.U.lLUJZ 
oOi>«ou*^aeuj 

■ 93 9 - 


O O o o o u 


U O 0-0 u o o 


ooouoouu 


<OKC 09 ‘ 0 >^'M(*^ 

OOOOOOOO 

oooooooo 

zzzzzzzz 


fom<ncn«n444 

oooooooo 

oooooooo 


04tf^'0^>coO‘0^(Mr<>4l^^4coO‘ 

44444't4iA'Aini/>U>i/>lf\in«/N 

oooooooooooooooo 

ooooooooooooooop 

zzzzzzzzzzzzzzzz 


85 


Figure 27 SOURCE AND CROSS REFERENCE LISTINGS FOR SUBROUTINE CONVLV 
(PAGE 3 OF 4) 
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expressed mathematically as 

= f ^ (i) f (28) 

the first or primary function in the calling list must contain more points 
than the second function. If this is not the case the program returns zeros 
for all points in the convolved function. Equation 28 is implemented as 


f- 


f = i- f 


where = 0.0 for » ^ o 

71 


or rt ^ rtc + 1 f 

^ c ' + t 


f- , - Q . o • The assumption that the primary function is zero outside of 
the rcinge of its input should be compensated for in the calling program if it 
is not valid. This is accomplished by deleting points at the beginning and/or 
the end of the convolved function array returned to the calling program. Note 
that the total number of points in the returned function is one more than the 
sum of the number of points in each individual input function. 


When either the differentiation or smoothing operations have been 
specified Equation 28 is implemented as 


=Z / 9r> J ' 

?» =/ t V- ri- 


» ' • ' 7 


^ '^c, 2 = • * ' = 

2 


c,v^ c,7.^ — 


where is set equal to the appropriate internally stored filter and 

is the number of points in the filter array. Note that points, 

divided equally between the beginning and the end of the convolved function 
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are not computed by Equation 30 but arbitrarily set to the nearest computable 
point. Therefore, although points are returned in the convolved 

function array, some of these points contain artificial information and may 
require special treatment in the calling program. 

j 

The major difference between Equations 29 and 30, beside the method 
of treatment of the points at the beginning and end of the convolved function 
array, is the folding of the secondary function ^ . In Equation 29 this 

folding is accomplished through manipulation of the index while this is not 
required in Equation 30 because of symmetry properties of the internally 
stored smoothing and differentiating filters. The secondary fimction is 
actually set equal to the negative of the stored differentiating filter in 
order to accomplish the required folding. 

Note that both Equation 29 and Equation 30 do not account for the 
effect of the differential, , present in -Equation 28. This must be 
accomplished in the calling program if required. 

5.0 Subroutine DXSCAL - This program is part of the Edge Gradient Spectral 

Analysis software package. It supports Program EGSA in the calculation 
of the OTF from edge trace data. The purpose of this subroutine is to determine 
the best scale factor for the x-axis for plots of individual functions 
requested by the user as an output option. No direct input data must be 
furnished by the user. The input data is supplied through the calling program. 
No printed output data is furnished by this program. 

No explicit listing of this program is provided since its require- 
ments may depend upon the user's computing facility similar to the other 
plotting subroutines. Its operation, however, is straightforward and will 
peimit a user to code an equivalent program. 

■Haree calls to DXSCAL are made by Program EGSA. The calling sequence 
is DXSCAL CXL, NP, DX, DXT, DXI) where 
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XL = the length of the x-axis in inches (input variable) 

NP = the number of points in the array to be plotted (input variable) 
DX = the increment between points (input variable) 

DXT = the total length of the x-axis in units of DX (output variable) 
DXl = the scale for plotting in inches/point (output variable) 

The purpose of the subrouting is to compute "convenient" values for DXT and 
DXl. DXl is chosen so that DX/DXI is the smallest number in the set (10™, 
1.5*10™, 2*10™, 2.5*10™, 4*10™, 5*10™, 8*10™, m an integer) that makes the 
data fit an axis of length XL, i.e. subject to the restriction that (NP-1) * 

DXl i XL. This scaling ensures that the plotted function will fill at least 
half of the XL inches. 

DXl is used directly by Program EGSA to compute an x-axis array for 
the plots of individual edges, line spreads or smoothed line spreads. The 
DXT variable is placed on the plots so that the axis can be labeled by user 
if desired. The units of DX and hence DXT in these cases are mm. 

6.0 Subroutine FORINV - This program is part of the Edge Gradient Spectral 

Analysis software package. It supports Program EGSA in the calculation of 
an OTP from edge trace data. The purpose of this subroutine is to Fourier 
transform complex spectra (modulus and phase) to obtain a real output 
function. Specifically it is used by Program EGSA to compute the smooth line 
spread function from the truncated OTF. No direct input data must be furnished 
by the user. The input data is supplied through the calling program. No 
printed output data is furnished by this program. 

Figure 28 is a listing of subroutine FORINV. The function of the 
major blocks are identified by the comments included in the coding. The 
algorithm employed in this program is specifically designed for transforming 
a complex to a real function. The function to be transformed is expressed 
as modulus and phase functions. Consequently the inverse Fourier transform 
can be written as 

I(x) = Jz(y)e clu (31a) 
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(31b) 


and the discrete representation as 

n = -n* 

where r^=r^v„), <p rt “ ^ ) » ^“id j^„ = - . Note that the effect of 

the differential, c/P , in Equation 31a has been suppressed in Equation 31b. 
Since I(x ) is real the following conditions must be met 

0.0 (32) . 

Expanding Equation 31b and using these conditions we can write the discrete 
form of the inverse transform as 


r M = Y. ^r, ( 2 n p^ 

^) + L Br, ^t^( 2 n p.„x) 

(33) 

n-0 

r>>^0 


where , Ay, ^ 

, 8 ^=>0.0 , By,=^ 2Z yy ^ „ 

and the 


factor of 2 accounts for the sununation over the negative values of n . For 
coding purposes this equation is rewritten with summation over positive 
indices by letting i - , namely , 

^ X ^Zir pf x) (34) 

i^i i=f 


and the coefficients Ai and are appropriately redefined. 

Xj ‘ C-^ i /=/;"•; 2 ^ 

^ 1 = - ')^y i ^ f f ' " > 77*+! 


Letting 

(35a) 

(35b) 


and requiring that 


we can write 



(35c) 
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(35d) 


21T Uj Xj 


TT(i-l) (-n*+j) 
rt * 


Setting . = (i-i) (- Ti*y) 

transform can be written^ 


in Equation 35d, 






I(xj) .1 "ai ^ (^) . £ 

t - / < = / 


7)*+ 1 


B ^ 




(36) 


The sine and cosine terms need only be computed over the first quadrant, i.e. 

y,* 

O i y» ii £ — ^ appropriate algorithm developed for selecting the 

a a 

appropriate element from this array when the argument lies in another quadrant 
Note that w/y can be less than zero and is accounted for in the coding by 
changing the sign of the second summation in Equation 36. In order to select 
the appropriate sine or cosine term and identify the quadrant the initial 
step involves a transformation of from the range [o, to [o , 2 r?'^) 

through the expression 

I = w.- • - 2 rt* 


m. 


I 2 77* 


(37) 


where I 1 represents the greatest integer which is less than the 
enclosed argument. The quadrant can then be subsequently specified by the 
expression 



o< ^ 


(38) 


and the value of L given by Equation 37 appropriately modified depending 
upon the quadrant to select the appropriate sine and cosine terms to perform 
the summations indicated by Equation 36. 


The restriction of performing an inverse transform from a complex 
to a real function, i.e. conditions of Equation 32 and the required relationship 
between the spatial and frequency increments. Equation 35c, permit the develop- 
ment of this algorithm which results in a considerable savings in computation 
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time. As the listing in Figure 28 indicates the user has an option of smoothing 
the modulus prior to inverse transformation. If this option is elected the 
initial step in the program is the setup of the smoothing factors SIGF(I). This 
is followed by the calculation of the coefficients A(I) and B(I) identified at 
Equation 33. The sine and cosine factors appearing in Equation 36 are then 
computed for argument values in the first quadrant and the calculation of both 
sum terms of the inverse transform of Equation 36 completed. In calculating 
the sums the program branches depending upon the value of KSIGN which is 
equivalent to Equation 38. After all the sum terms are computed, the spatial 
increment is defined using Equation 35c and execution is returned to the 
calling program. 

7.0 Subroutine FOURTR - This program is part of the Edge Gradient Spectral 

Analysis software package. It supports Program EGSA in the calculation of an 
OTF from edge trace data. The purpose of this subroutine is to Fourier trans- 
form a real input function to obtain a complex spectrum expressed by modulus 
and phase functions. It is specifically employed by Program EGSA to calculate 
individual OTF from differentiated edge trace data. No direct input data must 
be furnished by the user. The input data is supplied through the calling 
program. No printed output data is furnished by this program. 


Figure 29 is a listing of subroutine FOURTR. The function of the 
major blocks are identified by the comments included in the coding. The 
algorithm employed in this program is specifically designed for transforming 
a real to a complex function and results in considerable savings in execution 
time compared to a more direct implementation. If we assume that the trans- 
formed function will be written in terms of a modulus function T (u) and a 
phase function <(>(v) the Fourier transform can be expressed by the equation 

iS^u) 2tT iv X 

I(x)e dx C39a) 

-00 ■ 


The corresponding discrete representation 


t(u) 

-n* 


2TT i VXyj 


77= -T 



(39b) 



this subroutine is part of the edge gradient spectral analysis 
lEGSAI software package, this package estihates the optical 
transfer function' (otfi for any optical system using a set of edge 
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DELETE EXECUTION IF "NDATA" > 1000 
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Figure 29 SOURCE AND CROSS REFERENCE LISTINGS FOR SUBROUTINE FOURTR 
(PAGE 2 OF 4) 
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Figure 29 SOURCE AND CROSS REFERENCE LISTINGS FOR SUBROUTINE FOURTR 
(PAGE 3 OF 4) 
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where %^= 0.0 , increment between points. Note 

that the effect of the differential, , is included in the discrete repre- 
sentation. Since I (x) is real, the summation can be divided into real and 
imaginery terms, namely 


T: Cv) Q + i b(v) 


(40a) 


where 


and 




7) = - V 


CU7X} (2-rtV Xy,) 


(40b) 


q(v) = a^Y. Z TT Xyj ) 

n = - n* 


(40c) . 


Through algebraic manipulation the sum over the negative and positive indices 
can be combined and for coding purposes the summation rewritten using positive 
indices by letting j =^ 7 >+f ; j ^ f , • • • . Equations 40b and 40c become 


where 


•n^+f 

A (p) = F - (zTT F Xj) 

(41) 



■n^ + 1 


B(u)=Aj^Y. 

(42) 


/=' 



I(o)} 3 = i 


(43) 




0.0 ; y = / 

I(xj) - I(-Xj) ; j>z 


(44) 
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The modulus and phase of the Fourier transform can be written in terms of 
the real and imaginery terms by referring to Equation 40a, namely 


rM = 

(a ^ b(u)‘ ) 

-(45) 

(p(v) = 

^ A(t>) 

(46) 


The arguments of the sine and cosine factors, appearing in Equations 
41 and 42 are treated in the same fashion described in the previous section, 
namely setting 



(47a) 

i = f , • * ' ; 77^+ f 

(47b) 


and requiring that 




(47c) 


we can write the argument of the sine and cosine terras as 


2 7T Xj 


7T(i-l)(}-f) ^ 7T ynjj 


7 ) 


(47d) 


In this case, however, is not negative and has values in the range 

. The equations for the real and imaginery terras become 



The sine and cosine terms need only be computed over the first quadrant, i.e. 
O ^ mij ^ — , and an appropriate algorithm developed for selecting the 
correct element from this array when the argioraent lies in another quadrant. 
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We define 



(50) 


where I 1 : represents the greatest integer which is less than the 

enclosed argioment. The quadrant is then specified by the expression 

K = 1 "^] +1 i o ± K 6 ( 51 ) 

and the value of L from Equation 50 appropriately modified depending upon 
the quadrant to select the correct sine and cosine terms from the stored 
arrays . 


The restriction of Fourier transforming a real function to obtain 
a complex spectrum and the required relationship between spatial and frequency 
increments. Equation 47c, are required to develop the algorithm described and 
permit a considerable savings in computation time. The elements of the 
modulus and phase arrays are calculated using a modified form of Equations 45 
and 46, namely 


T. = (A^ + 


(f>- = 




(52) 

(53) 


Note that the phase array is equivalent to that described by Equation 46. 
However, the modulus has been multiplied by the spatial frequency increment 
and appropriate action should be taken in the calling program if 
required. 

The first step in the program is to adjust the number of points in 
the input function array to be even. This is not a requirement of the 
mathematics but executed only as an interface convenience with subroutine 
FOUINV described in the previous section. The origin of the spatial axis. 
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i.e. X-0 , is assumed to be located at the midpoint of the array. The next 
step in the program is the establishment of the sum and difference arrays 
expressed by Equations 43 and 44 and subsequently used to calculate the real 
and imaginary parts of the Fourier transform. This is followed by the 
computation of the sine and cosine teims within the first quadrant. The real 
and imaginery terms can then be computed from Equations 48 and 49 and subsequent 
branching depending upon the quadrant number specified by the variable KSIGN 
calculated using Equation 51. The sign of the sine and cosine terms are 
adjusted in accordance with the quadrant number and selected from the 
corresponding array using the transformed value of L (c.f. Equation 50). 

The final step in the program is the calculation of the modulus and phase 
arrays using Equations 52 and 53 and the definition of the spatial frequency 
increment based upon Equation 47c. Execution is then returned to the calling 
program. 

8.0 Main Program TONEQ - This program provides data which can be used to 
evaluate the tone quality of a photograph. It computes the gain achieved as 
a function of object contrast and average exposure level. Input data on 
punched cards must be supplied by the user. Printed values of the gain factors 
are furnished in tabular form by the program. 

Figure 30 is a listing of TONEQ. The function of the major blocks 
are identified by the comments included in the coding. In order to exercise 
this program two subprograms that are part of the Calspan Corporation Program 
Library are required. Detailed descriptions of these programs will not be 
furnished in this Appendix, however, the user is likely to have equivalent 
programs available in the library at his computing facility. CLEAR fills 
each byte of a specified block of storage with zeros. This operation can be 
accomplished by substituting a sequence of do loops in the TONEQ program. 

SPLN46 calculates a set of interpolation functions for a set of input points. 
These interpolation functions produce a continuous function that has a 
continuous first and second derivative over the interval spanned by the input 
data. The function tends to avoid waviness often encountered by other curve 
fitting techniques. The input data points must be reasonably accurate or 
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Figure 30 SOURCE AND CROSS REFERENCE LISTINGS FOR MAIN PROGRAM TONEO 
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OBTAIN SPLINE INTERPOLATION FUNCTIONS FOR SENSITONETRIC CALIBRATION 
DATA AND PRINT OUT COEFFICIENTS IN TABULAR FORM 
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SYMBOL 1NTEI«NAL STATEMENT NUMBERS 
U 0002 0003 0015- 0023 0025 0037 0042- 

I 00A4 OOA^ 00^ 00A5 0046 0047 
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Figure 30 SOURCE AND CROSS REFERENCE LISTINGS FOR MAIN PROGRAM TONEQ 
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smooth since the subroutine forces the function to pass through these points 
and consequently cannot perform any smoothing of the input data. The inter- 
polation functions are piecewise cubic and yield about the same curve one 
would obtain by forcing a flexible straight edge (a spline) through the set 
of input points. For additional information about spline interpolation 
procedures the reader is referred to DeBoor, Carl; Bicubic Spline Interpolation 
Journal of Mathematics and Physics , Vol. 41, No. 3, March 1962, pp. 212. 
Subroutine SPLN46 is used by this program to obtain interpolation functions 
for the sensitometric calibration data. 


The initial step in the program is to read user input data which 
specifies the number of points in the sensitometric calibration data, the 
sense of the imagery (i.e. positive or negative) and the number from 1 to 21 
of the first step in the calibration data. The user may also include a title 
to identify the photograph under evaluation if he so desires. Figure 31 shows 
a typical input data set. This example is taken from that used in the 
evaluation of the Apollo 17 panoramic photography. The first two cards compose 
the initial data block and the first card indicates that there are 14 points 
in the sensitometric calibration data measured on a negative image. The 
number of the step for the first data point is 8. The fact that the magnitude 
of the variable SENSE is 1.0 indicates that the difference in relative log 
exposure between points in the sensitometric calibration is equal to 0.15. 

The optional identification label indicates that the data were obtained from 
a Direct Negative roll of Apollo 17 panoramic photography that contains 
Frame 3111. An optional identification tag can be placed in columns 73 through 
80 to aid in data set identification. The figure shows a second data block for 
another photograph; in this case a Third Generation Positive copy of the same 
frame. This can be followed by additional data blocks for each photograph 
under evaluation. 

The next step in the program is the input of the sensitometric 
calibration value (density values) . This data is contained on the second 
card of the first data block and the second and third cards in the second data 
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Figure 31 TYPICAL TONEQ PROGRAM INPUT DATA SET 


block. The program then calculates an array of relative log exposure values 
X(N) in the range 0.0 to 3.0 that correspond to the input density values. 

By calling subroutine SPLN46 the coefficients of the spline interpolation 
functions for the sensitometric calibration data are obtained and subsequently 
printed out in tabular form. Figure 32 shows the printout generated for the 
first data block of Figure 31. The next step in the program involves the 
computation of the gain factors at an average exposure level equal to that of 
each point in the calibration data and three different object contrast levels. 
The levels of object contrast are specified by a data statement at the beginning 
of the program. Once the calculations have been completed the gain factors 
are printed out in tabular form. An example of this printout is included in 
Figure 32 for the first data block identified in Figure 31. The program then 
recycles to process the next data block if one is supplied by the user. If 
not, the program terminates execution. 
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Figure 32 EXAMPLE OF TONEQ PROGRAM PRINTED OUTPUT 
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